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Abstract: In these lectures, we develop the theory of the Colour Glass 
Condensate. This is the matter made of gluons in the high density environ- 
ment characteristic of deep inelastic scattering or hadron-hadron collisions 
at very high energy. The lectures are self contained and comprehensive. 
They start with a phenomenological introduction, develop the theory of 
classical gluon fields appropriate for the Colour Glass, and end with a 
derivation and discussion of the renormalization group equations which 
determine this effective theory. 



1 General Considerations 
1.1 Introduction 

The goal of these lectures is to convince you that the average properties of 
hadronic interactions at very high energies are controlled by a new form of mat- 
ter, a dense condensate of gluons. This is called the Colour Glass Condensate 
since 

• Colour: The gluons are coloured. 



Glass: The associated fields evolve very slowly relative to natural time 
scales, and are disordered. This is like a glass which is disordered and is 
a liquid on long time scales but seems to be a solid on short time scales. 



•Lectures given at the NATO Advanced Study Institute "QCD perspectives on hot and 
dense matter", August 6-18, 2001, in Cargese, Corsica, Prance 
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• Condensate: There is a very high density of massless gluons. These gluons 
can be packed until their phase space density is so high that interactions 
prevent more gluon occupation. With increasing energy, this forces the 
gluons to occupy higher momenta, so that the coupling becomes weak. 
The gluon density saturates at a value of order l/ag ^ 1, corresponding 
to a multiparticle state which is a Bose condensate. 

In these lectures, we will try to explain why the above is very plausible. 
Before doing this, however, it is useful to review some of the typical features 
of hadronic interactions, and some unanswered theoretical questions which are 
associate with these phenomena. This will motivate much of the later discussion. 

1.2 Total Cross Sections at Asymptotic Energy 

Computing total cross sections as i? ^ oo is one of the great unsolved problems 
of QCD. Unlike for processes which are computed in perturbation theory, it is 
not required that any energy transfer become large as the total collision energy 
E ^ oo. Computing a total cross section for hadronic scattering therefore 
appears to be intrinsically non-perturbative. In the 60's and early 70's, Regge 
theory was extensively developed in an attempt to understand the total cross 
section. The results of this analysis were to our mind inconclusive, and certainly 
can not be claimed to be a first principles understanding from QCD. 
The total cross section for pp and pp collisions is shown in Fig. ^ Typically, it 
is assumed that the total cross section grows as In^ E as E ^ oo. This is the 
so called Froissart bound, which corresponds to the maximal growth allowed 
by the unitarity of the scattering matrix. Is this correct? Is the coefficient of 
In^ E universal for all hadronic precesses? Why is the unitarity limit saturated? 
Can we understand the total cross section from first principles in QCD? Is it 
understandable in weakly coupled QCD, or is it an intrinsically non-perturbative 
phenomenon? 

1.3 Particle Production in High Energy Collisions 

In order to discuss particle production, it is useful to introduce some kinematical 
variables adapted for high energy collisions: the light cone coordinates. Let 
z be the longitudinal axis of the collision. For an arbitrary 4-vector = 
{v^ ,v'^ ,v^) {v^ = Vz, etc.), we define its light-cone (LC) coordinates as 

v+^^{v° + v'), v-^^iv'^-v'), v^^{v\v^). (1.1) 

In particular, we shall refer to = {t + as the LC "time", and to 

~ (i — z)/'\/2 as the LC "longitudinal coordinate" . The invariant dot product 
reads: 

p ■ X ^ p^x^ +p'^x^ — p± ■ x±, (1-2) 
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Figure 1: The cross sections for pp and pp scattering. 
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which suggests that p~ — the momentum variable conjugate to the "time" x"*" 
— should be interpreted as the LC energy, and p"*" as the (LC) longitudinal 
momentum. In particular, for particles on the mass- shell: _p± = {E±p^)/V2, 
with E = {m? + p^)^/^, and therefore: 

= \{E^-pI)^\{pI+'^^) = 1^1- (1-3) 
This equation defines the transverse mass m± . We shall also need the rapidity : 

y^h.E;.y±^. ,,4, 

2 p 2 TOj_ 

These definitions are useful, among other reasons, because of their simple prop- 
erties under longitudinal Lorentz boosts: p"*" —f np'^, p~ —^ {1/k)p~ , where k is 
a constant. Under boosts, the rapidity is just shifted by a constant: y ^ y-|- k. 
Consider now the collision of two identical hadrons in the center of mass frame, 
as shown in Fig. In this figure, we have assumed that the colliding hadrons 
have a transverse extent which is large compared to the size of the produced 
particles. This is true for nuclei, or if the typical transverse momenta of the 
produced particles is large compared to Aqcd, since the corresponding size will 
be much smaller than a Fermi. We have also assumed that the colliding particles 
have an energy which is large enough so that they pass through one another and 
produce mesons in their wake. This is known to happen experimentally: the 
particles which carry the quantum numbers of the colliding particles typically 
lose only some finite fraction of their momenta in the collision. Because of their 
large energy, the incoming hadrons propagate nearly at the speed of light, and 
therefore are Lorentz contracted in the longitudinal direction, as suggested by 
the figure. 

In LC coordinates, the right moving particle ("the projectile") has a 4- momentum 
Pi ^ {Pi ■,Pi :^±) with p^ ~ V2pz and p^ — A'P/2pf (since Pz 3> M, with 
M = the projectile mass). Similarly, for the left moving hadron ("the tar- 
get"), we have = Pi and P2 = Pi . The invariant energy squared is 
s = {pi + P2)^ ~ 2pi ■ p2 — 2piP2 — 4p^, and coincides, at it should, with 
the total energy squared {Ei + £^2)^ in the center of mass frame. 
We define the longitudinal momentum fraction, or Feynman's x, of a produced 
pion as 

X ^ 4 (1-5) 
Pi 

(with < X < 1). The rapidity of the pion is then 

y = 7:ln — = - In — — ^ Vproj - In - + In , (1.6) 

2 p^ 2 X TOj_ 

where jproj — ln[^/2p'^/M) = ln(Y^/Af). The pion rapidity is in the range 
—Jproj ^ y S5 Yproj (up to an overall shift by Ay — ln(Af/mj_)). 
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Figure 2: A hadron-hadron collision. The produced particles are shown as 
circles. 
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Figure 3: The rapidity distribution of particles produced in a hadronic collision. 



A typical distribution of produced particles (say, pions) in a hadronic collision 
is shown in Fig. ^. We denote by dN / dy the number of produced particles per 
unit rapidity. The leading particles are shown by the solid line and are clustered 
around the projectile and target rapidities. For example, in a heavy ion collision, 
this is where the nucleons would be. The dashed line is the distribution of 
produced mesons. Several theoretical issues arise in multiparticle production: 
Can we compute dN/dy ? Or even dN/dy at y = ("central rapidity") ? 
How does the average transverse momentum of produced particles {p±) behave 
with energy? What is the ratio of produced strange/nonstrange mesons, and 
corresponding ratios of charm, top, bottom etc at y = as the center of mass 
energy approaches infinity? Does multiparticle production as s — > cx) at y = 
become simple, understandable and computable? 

Note that y = corresponds to particles with = or p+ = m±/\/2, for 
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dN 

dy 




Figure 4: Feynman scaling of rapidity distributions. The two different lines 
correspond to rapidity distributions at different energies. 



which X = m^/ {\/2pi) = m^j^fs is small, x ^ 1, in the high-energy limit 
of interest. Thus, presumably, the multiparticle production at central rapidity 
reflects properties of the small-x degrees of freedom in the colliding hadron 
wavefunctions. 

There is a remarkable feature of rapidity distributions of produced hadrons, 
which we shall refer to as Feynman scaling. If we plot rapidity distributions 
of produced hadrons at different energies, then as function of y — Yproj i 

the 

rapidity distributions are to a good approximation independent of energy. This 
is illustrated in Fig. ^, where the rapidity distribution measured at one energy 
is shown with a solid line and the rapidity distribution at a different, higher, 
energy is shown with a dotted line. (In this plot, the rapidity distribution at 
the lower energy has been shifted by an amount so that particles of positive 
rapidity begin their distribution at the same yproj as the high energy particles, 
and correspondingly for the negative rapidity particles. This of course leads to 
a gap in the center for the low energy particles due to this mapping.) 
This means that as we go to higher and higher energies, the new physics is 
associated with the additional degrees of freedom at small rapidities in the center 
of mass frame (small-x degrees of freedom). The large x degrees of freedom do 
not change much. This suggests that there may be some sort of renormalization 
group description in rapidity where the degrees of freedom at larger x are held 
fixed as we go to smaller values of x. We shall see that in fact these large x 
degrees of freedom act as sources for the small x degrees of freedom, and the 
renormalization group is generated by integrating out degrees of freedom at 
relatively large x to generate these sources. 

1.4 Deep Inelastic Scattering 

In Fig. ||, deep inelastic scattering is shown. Here an electron emits a virtual 
photon which scatters from a quark in a hadron. The momentum and energy 
transfer of the electron is measured, but the results of the hadron break up are 
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Figure 5: Deep inelastic scattering of an electron on a hadron. 



not. In these lectures, we do not have sufficient time to develop the theory 
of deep inelastic scattering (see, e.g., for more details). For the present 
purposes, it is enough to say that, at large momentum transfer 3> ^cd^ 
this experiment can be used to measure the distributions of quarks in the hadron. 
To describe the quark distributions, it is convenient to work in a reference frame 
where the hadron has a large light-cone longitudinal momentum P+ 3> M 
("infinite momentum frame"). In this frame, one can describe the hadron as 
a collection of constituents ("partons"), which are nearly on-shell excitations 
carrying some fraction x of the total longitudinal momentum . Thus, the 
longitudinal momentum of a parton is p"*" = xP"*", with < x < 1. 
For the struck quark in Fig. ^, this x variable ("Feynman's x") is equal to the 
Bjorken variable x^j, which is defined in a frame independent way as x^j = 
/2P ■ q, and is directly measured in the experiment. In this definition, = 
—qf^Qf^, with the (space-like) 4-momentum of the exchanged photon. The 
condition that x = x^j is what maximizes the spatial overlap between the 
struck quark and the virtual photon, thus making the interaction favourable. 
The Bjorken variable scales like xbj ~ Q'^/s, with s = the invariant energy 
squared. Thus, in deep inelastic scattering at high energy (large s at fixed Q^) 
one measures quark distributions dN quark /dx. at small x (x <C 1). 
It is useful to think about these distributions as a function of rapidity. We define 
the rapidity in deep inelastic scattering as 

y = Yhadron - ln(l/x), (1.7) 

and the invariant rapidity distribution as 

dN dN 

In Fig. ^, a typical dN/dy distribution for constituent gluons of a hadron is 
shown. This plot is similar to the rapidity distribution of produced particles in 
hadron-hadron collisions (see Fig. ^). The main difference is that, now, we have 
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only half of the plot, corresponding to the right moving hadron in a collision in 
the center of mass frame. 

One may in fact argue that there is indeed a relationship between the structure 
functions as measured in deep inelastic scattering and the rapidity distributions 
for particle production. We expect, for instance, the gluon distribution function 
to be proportional to the pion rapidity distribution. This is what comes out in 
many models of particle production. It is further plausible, since the degrees of 
freedom of the gluons should not be lost, but rather converted into the degrees 
of freedom of the produced hadrons. 

The small x problem is that in experiments at HERA, the rapidity distributions 
for quarks and gluons grow rapidly as the rapidity difference 

T = ln(l/x) = Yhadron - J (1.9) 

between the quark and the hadron increases [||. This growth appears to be 
more rapid than t or r^, and various theoretical models based on the original 
considerations by Lipatov and colleagues ^ suggest it may grow as an expo- 
nential in r H, The more established DGLAP evolution equation predicts 
a less rapide growth, like an exponential in -y/r, but this is still exceeding the 
Froissart unitarity bound, which requires rapidity distributions to grow at most 
as (since r '--^ Ins). 

In Fig. 0, the ZEUS data for the gluon distribution are plotted for = 
5 GeV^ 20 GeV^ and 200 GeV^ ^. The gluon distribution is the number of 
gluons per unit rapidity in the hadron wavefunction, xG'(x, Q^) — dNgiuons/dy . 
Experimentally, it is extracted from the data for the quark structure functions, 
by exploiting the dependence of the latter upon the resolution of the probe, that 
is, upon the transferred momentum . Note the rise of xG(x, Q^) at small x: 
this is the small x problem. If one had plotted the total multiplicity of produced 
particles in pp and pp collisions on the same plot, one would have found rough 
agreement in the shape of the curves. 
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Figure 7: The Zeus data for the gluon structure functions. 




Figure 8: Saturation of gluons in a hadron. A view of a hadron head on as x 
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Why is the small x rise in the gluon distribution a problem? Consider Fig. ||, 
where we view the hadron head on. The constituents are the valence quarks, 
gluons and sea quarks shown as coloured circles. As we add more and more 
constituents, the hadron becomes more and more crowded. If we were to try 
to measure these constituents with say an elementary photon probe, as we do 
in deep inelastic scattering, we might expect that the hadron would become 
so crowded that we could not ignore the shadowing effects of constituents as 
we make the measurement. (Shadowing means that some of the partons are 
obscured by virtue of having another parton in front of them. This would result 
in a decrease of the scattering cross section relative to what is expected from 
incoherent independent scattering.) 

We shall later argue that the distribution functions at fixed saturate, which 
means that they cease growing so rapidly at high energy 0, |[ ^, |l^ . (See also 
Refs. [0 |l^, for recent reviews and more references.) This saturation 
will be seen to occur at transverse momenta below some intrinsic scale, the 
"saturation scale" , which is estimated as: 

Ql^a.N,-^^, (1.10) 

where dN/dy is the gluon distribution. Only gluons matter since, at small x, 
the gluon density grows faster then the quark density, and is the driving force 
towards saturation. This is why in the forthcoming considerations we shall 
ignore the (sea) quarks, but focus on the gluons alone. Furthermore, 7ri?^ — 
with R the hadron radius — is the area of the hadron in the transverse plane. 
(This is well defined as long as the wavelengths of the external probes are small 
compared to R.) Finally, UgNc is the colour charge squared of a single gluon. 
Thus, the "saturation scale" ( 1.1C| ) has the meaning of the average colour charge 



squared of the gluons in the hadron wavefunction per unit transverse area. 
Since the gluon distribution increases rapidly with the energy, as shown by the 
HERA data, so does the saturation scale. We shall use the rapidity difference 



T = ln(l/x) '-^ Ins, eq. (l.E), to characterize this increase, and write = Qgir) 



For sufficiently large r (i.e., high enough energy, or small enough x), 

Qlir) » AlcD, (1-11) 

and as{Q^) <C 1. Then we are dealing with weakly coupled QCD, so we should 
be able to perform a first principle calculation of, e.g., 

• the gluon distribution function; 

• the quark and heavy quark distribution functions; 

• the intrinsic p± distributions of quarks and gluons. 

But weak coupling does not necessarily mean that the physics is perturbative. 
There are many examples of nonperturbative phenomena at weak coupling. An 
example is instantons in electroweak theory, which lead to the violation of baryon 
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number. Another example is the atomic physics of highly charged nuclei, where 
the electron propagates in the background of a strong nuclear Coulomb field. 
Also, at very high temperature, QCD becomes a weakly coupled quark-gluon 
plasma, but it exhibits nonperturbative phenomena on large distances r 3> 1 /T 
(with T the temperature), due to the collective behaviour of many quanta pSf . 
Returning to our small-x gluons, we notice that, at low transverse momenta 
< Q1{t), they make a high density system, in which the interaction proba- 
bility 



a,(Q^) 1 dN 
Q2 7ri?2 dy 



1 when Q2 - gf (r) (1.12) 



is of order one 0, ^6). That is, although the couphng is small, as{Q^) <C 1, 
the effects of the interactions are amplified by the large gluon density (we shall 
see that dN / dy ^ 1 /as at saturation) , and ordinary perturbation theory breaks 
down. 

To cope with this, a resummation of the high density effects is necessary. Our 
strategy to do so — to be described at length in these lectures — will be to 
construct an effective theory in which the small-x gluons are described as the 
classical colour fields radiated by "colour sources" at higher rapidity. Physically, 
these sources are the "fast" partons, i.e., the hadron constituents with larger 
longitudinal momenta p"*" ^ xP"*" . The properties of the colour sources will be 
obtained via a renormalization group analysis, in which the "fast" partons are 
integrated out in steps of rapidity and in the background of the classical field 
generated at the previous steps. 

The advantage of this strategy is that the non-linear effects are dealt with in a 
classical context, which makes exact calculations possible. Specifically, (a) the 
classical field problem will be solved exactly, and (b) at each step in the renor- 
malization group analysis, the non-linear effects associated with the classical 
fields will be treated exactly. On the other hand, the mutual interactions of the 
fast partons will be treated in perturbation theory, in a "leading-logarithmic" 
approximation which resums the most important quantum corrections at high 
energy (namely, those which are enhanced by the large logarithm ln(l/x)). 
As we shall see, the resulting effective theory describes the saturated gluons as 
a Colour Glass Condensate. The classical field approximation is appropriate for 
these saturated gluons, because of the large occupation number Nk ^ I /as S> 1 
of their true quantum state. In this limit, the Heisenberg commutators between 
particle creation and annihilation operators become negligible: 

[flfc, 4] = 1 < alok = Nk , (1.13) 

which corresponds indeed to a classical regime. The classical field language 
is also well adapted to describe the coherence of these small-x gluons, which 
overlap with each other because of their large longitudinal wavelengths. 
The phenomenon of saturation provides also a natural solution to the unitarity 
problem alluded to before. We shall see that, with increasing energy, the new 
partons are produced preponderently at momenta p± Qs- Thus, these new 
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partons have a typical transverse size ^ l/p± ^ 1/Qs- Smaller is x (i.e., larger 
is t), larger is Qs(t), and therefore smaller are the newly produced partons. 
An external probe of transverse resolution Ax± ~ 1/Q will not see partons 
smaller than this resolution size. For r large enough, < Q'^{t), so that the 
partons produced when further increasing the energy will not contribute to the 
cross section at fixed Q^. Thus, although the gluon distribution keeps increasing 
with T, there is nevertheless no contradiction with unitarity. 



1.5 Geometrical Scaling 

Another striking feature of the experimental data at HERA is geometrical scal- 
ing at Bjorken x < 0.01 [|l^. In general, one expects the structure functions 
extracted from deep inelastic scattering to depend upon two dimensionless kine- 
matical variables, x and Q^/A^, where A^ is some arbitrary momentum scale 
of reference, which is fixed. The striking feature alluded to before is the obser- 
vation that the x dependence measured at HERA at x < 0.01 and for a broad 
region of g2 (between 0.045 and 450 Gev^) can be entirely accounted for by a 
corresponding dependence of the reference scale A^ 1/R^{x) alone. That is, 
rather than being functions of two independent variables x and Q^/A^, the mea- 
sured structure functions at x < 0.01 depend effectively only upon the scaling 
variable 

r = Q^R\x) (1.14) 

where (x) ~ x^ and A ~ 0.3 — 0.4 in order to fit the data. This is illustrated in 
Fig. [n^. Such a scaling behaviour is consistent with the saturation scenario 
[|8[ [l^, ]l^, as we shall discuss towards the end of these lectures. Note however 
that the experimentally observed scaling extends to relatively large values of 
X and Q^, above all the estimates for the saturation scale. Thus, this feature 
seems to be more general than the phenomenon of saturation. 



1.6 Universality 

There are two separate formulations of universality which are important in un- 
derstanding small X physics. 

a) The first is a weak universality j|, |l^. This is the statement that at suffi- 
ciently high energy, physics should depend upon the specific properties of the 
hadron at hand (like its size or atomic number A) only via the saturation scale 
Qs{t,A). Thus, at high energy, there should be some equivalence between nu- 
clei and protons: When their values are the same, their properties must 
be the same. An empirical parameterization of the gluon structure function in 
eq. ( p^ ) is 



1 dN A^l'^ 

(1.15) 



where (5 ^ 0.2 — 0.3 This suggests the following correspondences: 
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Figure 9: Experimental data on the cross section for virtual photon-proton deep 
inelastic scattering from the region x < 0.01 plotted verses the scaling variable 
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• RHIC with nuclei ~ HERA with protons; 

• LHC with nuclei ~ HERA with nuclei. 

Estimates of the saturation scale for nuclei at RHIC energies give ~ 1 — 2 Gev, 
and at LHC - 2 - 3 Gev. 

b) The second is a strong universality which is meant in a statistical mechanical 
sense. This is the statement that the effective action which describes small x 
distribution function is critical and at a fixed point of some renormalization 
group. This means that the behavior of correlation functions is given by uni- 
versal critical exponents, which depend only on general properties of the theory 
such as its symmetries and dimensionality. 



1.7 Some applications 

We conclude these introductory considerations with a (non-exhaustive) enumer- 
ation of recent applications of the concept of saturation and the Colour Glass 
Condensate (CGC) to phenomenology. 

Consider deep inelastic scattering first. It has been shown in Refs. that 
the HERA data for (both inclusive and diffractive) structure functions can be 
well accounted for by a phenomenological model which incorporates saturation. 
The same model has motivated the search for geometrical scaling in the data, 



as explained in Sect. 1.5 



Coming to ultrarelativistic heavy ion collisions, as experimentally realized at 
RHIC and, in perspective, at LHC, we note that the CGC should be the ap- 
propriate description of the initial conditions. Indeed, most of the multiparticle 
production at central rapidities is from the small-x (x < lO"*^) partons in the 
nuclear wavefunctions, which are in a high-density, semi-classical, regime. The 
early stages of a nuclear collision, up to times ^ l/Qs, can thus be described as 
the melting of the Colour Glass Condensates in the two nuclei. In Refs. , this 
melting has been systematically studied, and the multiparticle production com- 
puted, via numerical simulations of the classical effective theory ||, After 
they form, the particles scatter with each other, and their subsequent evolution 
can be described by transport theory [p2[ . 

The first experimental data at RHIC | |23{ have been analyzed from the perspec- 
tive of the CGC in Refs. Specifically, the multiparticle production 
has been studied with respect to its dependence upon centrality ("number of 
participants" ) , rapidity |^ and transverse momentum distribution Q . 
The charm production from the CGC in peripheral heavy-ion collisions has been 
investigated in |2^ . 

Electron-nucleus (eA) deeply inelastic scattering has been recently summarized 
in . Some implications of the Colour Glass Condensate for the central region 
of p + A collisions have been explored in Refs. ||2^, |3C[| . 
Instantons in the saturation environment have been considered in Ref. |3l|. 
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2 The classical effective theory 



With this section, we start the study of an effective theory for the small x 
component of the hadron wavefunction 
the previous review papers 
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(see also 



_ B B 

1). Motivated by the physical arguments 
exposed before, in particular, by the separation of scales between fast partons 
and soft (i.e., small-x) gluons, in the infinite momentum frame, this effective 
theory admits a rigourous derivation from QCD, to be described in Sect. ^j. 
Here, we shall rather rely on simple kinematical considerations to motivate its 
general structure. 



2.1 A stochastic Yang-Mills theory 

In brief, the effective theory is a classical Yang-Mills theory with a random 
colour source which has only a "plus" component Q : 

(D^F'^'^Ux) = 6^^+Pa{x). (2.1) 

The classical gauge fields represent the soft gluons in the hadron wavefunc- 
tion, i.e., the gluons with small longitudinal momenta (fc+ = xP+ with x <^ 1). 
For these gluons, the classical approximation should be appropriate since they 
are in a multiparticle state with large occupation numbers. 
The fast partons, with momenta p+ 3> fc"*", are not dynamical fields anylonger, 
but they have been rather replaced by the colour current = S'^'^pa which acts 
as a source for the soft gluon fields. This is quite intuitive: the soft gluons in 
the hadron wavefunction are radiated by typically fast partons, via the parton 
cascades shown in Fig. ^ It is in fact well known that, for the tree-level 
radiative process shown in Fig. ^a, classical and quantum calculations give 
identical results in the limit where the emitted gluon is soft 0|. What is less 
obvious, but will be demonstrated by the analysis in Sect. ^, is that quantum 
corrections like those displayed in Fig. |l^.b do not invalidate this classical 
description, but simply renormalize the properties of the classical source, in 
particular, its correlations. 

The gross properties of this source follow from kinematics. The fast partons 
move along the z axis at nearly the speed of light. They can emit, or absorb, 
soft gluons, but in a first approximation they preserve straightline trajectories 
along the light-cone (z = t). In terms of LC coordinates, they propagate in 
the positive direction, while sitting at = 0. Their colour current is 
proportional to their velocity, which implies = 5^^pa, with a charge density 
Pa{x) which is localized near x~ = 0. More precisely, as quantum fields, the 
fast partons are truly delocalized over a longitudinal distance Ax^ ^ 1/p^, as 
required by the uncertainty principle. But since <C 1/fc^, they still look as 
sharply localized when "seen" by the soft gluons, which have long wavelengths 
and therefore a poor longitudinal resolution. 

^Written as it stands, eg. ( |2.l| ) is correct only for field configurations hamn ct A~ = 0; when 
A~ ^ 0, the source p in its r.h.s. gets rotated by Wilson lines built from A~ |37|. 
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Figure 10: a) Soft gluon emission by a fast parton; b) a gluon cascade. 



The separation of scales in longitudinal momenta implies a corresponding sep- 
aration in time: Softer partons have larger energies, and therefore shorter life- 
times. Consider indeed the radiative process in Fig. |l^.a, where fc+ ^ This 
is a virtual excitation whose lifetime (in units of LC time x~^) can be estimated 
from the uncertainty principle as 

Ax+ = ~ — < — . (2.2) 

£p-k + £k — £p £k £p 

This is small as compared to the typical time scale 1 /sp for the dynamics of the 
fast partons. [In eq. ( ^.2| ), Ep = p^/2p+ is the LC energy of the on-shcU gluon 
with momentum p — {p^,pi_), and we have used the fact that, for <C p'^ 
and comparable transverse momenta k± and p±, £k ^ Sp,ep-k-] Thus, the 
"fast" degrees of freedom are effectively frozen over the short lifetime of the soft 
gluon, and can be described by a time-independent (i.e., independent of x~^) 
colour source Pa{x~ , x±). 

Still, this colour source is eventually changing over the larger time scale 1/Sp. 
Thus, if another soft gluon is emitted after a time interval ^ 1/ep, it will "see" a 
different configuration of p, without quantum interference between the different 
configurations. This can be any of the configurations allowed by the dynamics of 
the fast partons. We are thus led to treat pa{x^ , x±) as a classical random vari- 
able (here, a field variable), with some probability density, or weight function, 

[p\i which is a functional of p. 
As suggested by its notation, the weight function depends upon the soft scale 
fc"*" at which we measure correlations. Indeed, as we shall see in Sect. ^, Wk+ [p] 
is obtained by integrating out degrees of freedom with longitudinal momenta 
larger than k'^ . It turns out that it is more convenient to use the rapidit'^ 

r = ln(P+/fc+) =ln(l/x) (2.3) 

■^Strictly speaking, this is the rapidity difference between the small-x gluon and the hadron, 
as defined previously in eq. (IL^). But this difference is the relevant quantity for what follows, 
so from now on it will be simply referred to as "the rapidity" , for brevity. 
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to indicate this dependence, and thus write Wr[p] = Wk+[p]. 
To deal with field variables and functionals of them, it is convenient to consider 
a discretized (or lattice) version of the 3-diniensional configuration space, with 
lattice points {x~,x±). (We use the same notations for discrete and continu- 
ous coordinates, to avoid a proliferation of symbols.) A configuration of the 
colour source is specified by giving its values p°'{x' ,x±) at the N lattice points. 
The functional Wt-[/9] is a (real) function of these A'' values. To have a mean- 
ingful probabilistic interpretation, this function must be positive semi-definite 
(WV[p] > for any p), and normalized to unity: 



J D[p] Wr[p] = 1: 



with the following functional measure: 



(2.4) 



(2.5) 



Gluon correlation functions at the soft scale fc+ = xP+ — P+e are obtained 



by first solving the classical equations of motion {2A) and then averaging the 
solution over p with the weight function Wr[p\ (below x = {x^ ,xj_)) : 



{Al{x+,x)Ai{x+,y)- 



D[p] WAp]A:ix)Ai{y)^ 



(2.6) 



where Al^ = A\[p] is the solution to the classical Yang-Mills equations with 
static source Pa, and is itself independent of time (cf. Sect. 2.3 below). Note 
that only equal-time correlators can be computed in this way; but these are 
precisely the correlators that are measured by a small-x external probe, which 
is absorbed almost instantaneously by the hadron (cf. eq. ( |2.2| )). 
The formula (^^) is readily extended to any operator which can be related to p. 
To guarantee that only the physical, gauge-invariant, operators acquire a non- 
vanishing expectation value, we shall require WV[p] to be gauge-invariant. In 
practical calculations, one generally has to fix a gauge, so the gauge symmetry 
of Wt-[/o] may not be always manifest. 



To summarize, the effective theory is defined by eqs. {2A) and (2^) together 
with the (so far, unspecified) weight function Vl/V[p]- In what follows, we shall 
devote much effort to derive this theory from QCD, and construct the weight 
function Wr[p] in the process (in Sects. |^-J|). But before doing that, let us gain 
more experience with the classical theory by solving the equations of motion 
( ^?l| ) (in Sect. ^^), and th en u sing the result to compute the gluon distribution 
of a large nucleus (in Sect. 2A). In performing these calculations, we shall need 
a more precise definition of the gluon distribution function and, more generally, 
of the relevant physical observables, so we start by discussing that. 



2.2 Some useful observables 

In subsequent applications of the effective theory, we shall mainly focus on 
two observables which, because of their physical content and of the specific 
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structure of the effective theory, are particularly suggestive for studies of non- 
linear phenomena like saturation. These observables, that we introduce now, 
are the gluon distribution function and the cross-section for the scattering of a 
"colour dipole" off the hadron. 



2.2.1 The gluon distribution function 

We denote by G(x,Q^)d^ the number of gluons in the hadron wavcfunction 
having longitudinal momenta between xP^ and (x -|- (ix)P+, and a transverse 
size Ax± ~ 1/Q ■ In other terms, the gluon distribution xG{x, Q^) is the number 
of gluons with transverse momenta k± ^ Q per unit rapidity : 



xG{x,Q^) = J d^kj_k 



dN 



k+=xP+ 



dk+d^kl_ 

J d^k e(g2 - fci) x,5(x - k+/p+) ^ , (2.7) 



where k — (A;+,k_L) and 



dN dN 



d^k dk+d'^k. 



(2.8) 



is the Fock space gluon density, i.e., the number of gluons per unit of volume in 
momentum space. The difficulty is, however, that this number depends upon 
the gauge, so in general it is not a physical observables. Still, as we shall shortly 
argue, this quantity can be given a gauge-invariant meaning when computed in 
the light-cone (LC) gauge 

^+ = 0. (2.9) 

(We define the light-cone components of in the standard way, as = 
{Al ± A^)/y/2.) In this gauge, the equations of motion^ 

D^pt^'' = 0, (2.10) 

imply for the -|- component 

DiF'+ + D+P-+ ^0, (2.11) 
which allows one to compute A^ in terms of A^ as 

A- ^ ^D'd+A\ (2.12) 

This equation says that we can express the longitudinal field in terms of the 
transverse degrees of freedom which are specified by the transverse fields entirely 

•^For the purposes of LC quantization we use the equ atio ns of motion without sources; that 
is, we consider real QCD, and not the effective theory (2.1). 
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and explicitly. These degrees of freedom correspond to the two polarization 
states of the gluons. The quantization of these degrees of freedom proceeds by 
writing p9{ : 

J^,{x+.x) = j^^^^ (e*X(^+,^) (2.13) 

{x ■ k = x^k'^ — xj_ • k_L) with the creation and annihilation operators satisfying 
the following commutation relation at equal LC time a;+ : 

[a',{x+,k),a^^\x+,q)] = 5'^bc'2k+ {2^f S^'^Xk - q). (2.14) 

In terms of these Fock space operators, the gluon density is computed as: 

dN ^ ^ 9k+ 

^ = (af (x+, k) al{x+, k)) = ^ (A^(x+, fc)A^(x+, -k)) , (2.15) 

where the average is over the hadron wavefunction. By homogeneity in time, this 
equal-time average is independent of the coordinate x'^ , which will be therefore 
omitted in what follows. By inserting this into eq. ( |2.7| ) and using the fact that, 
in the LC-gauge, F^+(A:) = ik^ A\{k)^ one obtains (with fc+ = xP+): 

xG(x,Q2) = i j ^Q{Q^ -kl){Fl+{k)Fl+{-k)). (2.16) 

As anticipated, this docs not look gauge invariant. In coordinate space: 

F:+{k)Fl+{-k) = fd^x f d^ye^^^-y^-^^F:+{x;)F:+{y) (2.17) 



involves the electric fieldsQ at different spatial points x and y. A manifestly 
gauge invariant operator can be constructed by appropriately inserting Wilson 
lines. Specifically, in some arbitrary gauge, we define 

0-,{x,y) = Tr {F'+ {x) U^{x,y) F'+ {y) U^{y,x)] , (2.18) 
where (with Aa = (^+, A" ), Iee l„r°) 

C/,K,1.Pcxp{../..^l(,l}, (2,19) 

and 7 is an arbitrary oriented path from y to x. The (omitted) temporal coor- 



dinates x'^ are the same for all fields. For any path 7, the operator in cq. (2.18) 

is gauge-invariant, since the chain of operators there makes a closed loop. 

We now show that, by appropriately chosing the path, the gauge, and the bound- 



ary conditions, the gauge-invariant operator (2.18) can be made to coincide with 



*The component Fl'^ = — is usually referred to as the (LC) "electric field" by analogy 
with the standard electric field = F^^ = —d^A'^ (in the temporal gauge = 0). 



19 



A 



(- x^) 



(x, x^) 



(- ^, 



Figure 11: The path 7 used for the evaluation of the gauge- invariant operator 



the simple 2-point function ( 2.17 ). Specifically, consider the path shown in Fig. 
|lT| , with the the following three elements: two "horizontal" pieces going along 
the axis from {y^,yi_) to {—oo,y±), and, respectively, from {—oo,xj_) to 
{x~,x±), and a "vertical" piece from (— oo,j/_l) to {—co,x±). Along the hori- 
zontal pieces, dz- A = dz~ , so these pieces do not matter in the LC gauge. 
Along the vertical piece, dz- A = dz,± ■ Aj_{—oo, z±), and the path 7 between y± 
and x± is still arbitrary. But the contribution of any such a path to the Wilson 
line vanishes once we impose the following, "retarded" , boundary condition: 



as X 



(2.20) 



(Note that the "retardation" property refers here to x~ , and not to time.) 
To summarize, for the particular class of paths mentioned above, in the LC 
gauge ~ 0, and with the boundary condition (2.20), Uj{x,y) 1, and the 
manifestly gauge-invariant operator in eq. (2.18) reduces to the simpler operator 
( ^.17 ) which defines the number of gluons in this gauge. Converserly, the latter 
quantity has a gauge-invariant meaning, as the expression of a gauge-invariant 
operator in a specific gauge. 

We shall need later also the gluon distribution function in the transverse phase- 
space (in short, the "gluon density"), i.e., the number of gluons per unit rapidity 
per unit transverse momentum per unit transverse area: 



J\fr{k^,b^) = 



d^N 



dTd^kj_d%A 



d^ xG(x,fc2^) 

d^kj_d^bj_ 



(2.21) 



where r = ln(l/x) = ln(P^//j+) and b± is the impact parameter in the trans- 
verse plane (i.e., the central coordinate b± — {x± + y±)/2 in eq. (2.17)). This 



phase-space distribution is a meaningful quantity since the typical transverse 
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momenta we consider are relatively large, 



kl » ^QCD - ^/R\ (2.22) 

so that the transverse de Broglie wavelengths ~ of the partons under 

consideration are much shorter than the typical scale of transverse variation 
in the hadron, 1/Kqcd- (In particular, this explains why we can consider the 
hadron to have a well defined transverse size R.) 

In fact, for simplicity, we shall mostly consider a hadron which is homogeneous 
in the transverse plane, with a sharp boundary at radial distance R. Then, the 



eq. 



den sity (|2.2lD is independent of b±_ (within the disk b±_ < i?), and reads (cf. 



= ^ ^ = 4^ (fc)^r(-fc)). (2.23) 
2.2.2 The dipole-hadron cross-section 

Consider high energy deep inelastic scattering (DIS) in a special frame — the 
"dipole frame" — in which the virtual photon 7* is moving very fast, say, in the 
negative z direction, but most of the total energy is still carried by the hadron, 
which moves nearly at the speed of light in the positive z direction. Thus, the 
rapidity gap between the hadron and the virtual photon is 

T = yhadro7i - y-y- , with Ij-y'l Yhadron- (2.24) 



(As in Sect. 1.4, r ~ ln(l/x) sa \n{s/Q'^), where is the virtuality of 7* and 



s is the invariant energy squared. Note also that y^» < 0, since 7* is a left 
mover.) 

The dipole frame is special in two respects (and references therein): 

i) The DIS looks like a two step process, in which 7* fluctuates first into a 
quark-antiquark pair, which then scatters off the hadron. The qq pair is in a 
colour singlet state, so it forms a colour dipole. 

ii) The essential of the quantum evolution is put in the hadron wavefunction, 
which carries most of the energy. The dipole wavefunction, on the other hand, 
is simple and given by lowest order perturbation theory. More precisely, if 
as|y^.| <C 1, then the dipole is just a quark-antiquark pair, without additional 
gluons. 

Thus, in this frame, all the non-trivial dynamics is in the dipole-hadron scat- 
tering. Because of the high energy of the qq pair, this scattering can be treated 
in the eikonal approximation ^ ^| : the quark (and the antiquark) 

follows a straight line trajectory with z ~ —t (or = 0), and the effect of 
its interactions with the colour field of the hadronic target is contained in the 
Wilson line: 

V\x±) = Pexp(^ig J dx-A+ix-^x^)^^ , (2.25) 
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where x±^ is the transverse coordinate of the quark, t"'s are the generators of the 
colour group in the fundamental representation, and the symbol P denotes the 
ordering of the colour matrices A^ix) = A'^(x)t°' in the exponent from right to 
left in increasing order of their x~ arguments. Note that is the projection 
of along the trajectory of the fermion. For an antiquark with transverse 
coordinate y±_ the corresponding gauge factor is V{y±). Clearly, we adopt here 
a gau ge w here ^ (e.g., the covariant gauge to be discussed at length in 
Sect. 

It can then be shown that the S'-matrix element for the dipolc-hadron scattering 
is obtained by averaging the total gauge factor tviV'^ {x±_)V{y±_)) (the colour 
trace occurs since we consider a colourless qq state) over all the colour field 
configurations in the hadron wavefunction: 

Sr{x^,yi_) ^ ^(tr(0(x^)F(y^)))^. (2.26) 

The dipol e fram e is like the hadron i nfinit e momentum frame in that yhadron ~ 
T, cf . eq. ( 2.24 ) , so the average in eq. ( ^.26|) can be computed within the effective 



theory of Sect. 2A, that is, like in eq. (^.6|) 
The dipole-hadron cross section for a dipole of size r± = xj_ — y± is obtained by 
integrating 2(1 — Sr{xj_, yi_)) over all the impact parameters b± — {xj_ +yj_)/2: 



<Jd^pole{T,r^) = 2 J d^b^^(^tr(l-V^ix±)V{y±)'^)^. 



(2.27) 



Finally, the 7* -hadron cross-section is obtained by convoluting the dipole cross- 
section (2.27) with the probability that the incoming photon splits into a qq pair: 



<J^'h{T,Q^) ^ dzj d^r^\<^>{z,r^,Q^)\^'ad^pole{T,rJ^). (2.28) 

Here, 'i{z,r±;Q^) is the light-cone wavefunction for a photon splitting into a 
qq pair with transverse size r± and a fraction z of the photon's longitudinal 
momentum carried by the quark ^] . 

2.3 The classical colour field 

From the point of view of the effective theory, the high density regime at small 
X is characterized by strong classical colour fields, whose non-linear dynamics 
must be treated exactly. Indeed, we shall soon discover that, at saturation. 



xG'(x, Q^) ~ 1/q!s, which via eqs. (2.16) and (|2.6|) implies classical fields with 



amplitudes ~ 1/g. Such strong fields cannot be expanded out from the 
covariant derivative = 9* — ig A^. Thus, we need the exact solution to the 
classical equations of motion ( p.l| ), that we shall now construct. 
We note first that, for a large class of gauges, it is consistent to look for solutions 
having the following properties: 

f:'=0, A-=0, At,Ai: static, (2.29) 
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where "static" means independent of a;"*". (In fact, once such a static solution is 
found in a given gauge, then the properties ( 2.29| ) will be preserved by any time- 



independent gauge transformation.) This follows from the specific structure of 
the colour source which has just a "+" component, and is static. For instance. 



the component /i = i of eq. (2.1) reads 



= D^F""' = DjF^' + + D^F~\ (2.30) 



But_D+ = D = d — igA vanishes by eq. ( 2.29 ), and so does F Thu s 



eq. ( 2.3C ) reduces to DjF^^ — 0, which implies F*-' ~ 0, as indicated in eq. ( |2.2£ 
This further implies that the transverse fields A' form a two-dimensional pure 
gauge. That is, there exists a gauge rotation U{x^ ,x±) € SU(A^) such that (in 
matrix notations appropriate for the adjoint representation: A^ = A\T°-, etc) : 

A\x-,Xi_) = -U{x-,Xi_)d'U\x-,Xi_). (2.31) 
9 

Thus, the requirements ( |2.29| ) leave just two independent field degrees of free- 
dom, A'^{x) and U{x), which are further reduced to one (either yl+ or U) by 
imposing a gauge-fixing condition. 



We consider first the covariant gauge (COV-gauge) d^A>^ = 0. By eqs. (2.29) 



and (2.31), this implies diA^ = 0, or [/ = 0. Thus, in this gauge: 

i^(x) = 5^+aa{x-,xi_), (2.32) 

with aa{x) linearly related to the colour source pa in the COV-gauge : 

- Viaa(x) = ~pa{x). (2.33) 

Note that we use curly letters to denote solutions to the classical field equations 
(as we did already in eq. (|2.6|)). Besides, we generally use a tilde to indicate 
quantities in the COV-gauge, although we keep the simple notation aa{x) for 
the classical field in this gauge, since this quantity will be frequently used. 
Eq. ( ^.33 ) has the solution : 



aa{x ,xi_) = / d^2/_L (a;_L| — ^ |?/_L)pa(a; ,yi_) 

-r^li^T \2 2 Pa{x~,yi_), (2.34) 

47r (x±-y±yp^ 

where the infrared cutoff /i is necessary to invert the Laplacean operator in 
two dimensions, but it will eventually disappear from (or get replaced by the 
confinement scale Aqcd in) our subsequent formulae. 
The only non-trivial field strength is the electric field: 

= -d'aa. (2.35) 

In terms of the usual electric (E) and magnetic (B) fields, this solution is char- 
acterized by purely transverse fields, Ej^ = {E^,E^) and Bj^ = {B^,B^), which 
are orthogonal to each other: E_l • B_l = (since = —E"^ and B^ = E^). 
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To compute the gluon distribution ( 2.16| ), one needs the classical solution in 
the LC-gauge A+ = 0. Thi s is of the form Ai^ = with J^^{x- ,xa_) a 

The gauge rotation U{x) can be obtained by 
- to deduce an equation for U . 



"pure gauge", cf. eq. (2.31 



inserting the Ansatz (2.31) in eq. (2.1) with /i 



Alternatively, and simpler, the LC-gauge solution can be obtained by a gauge 



rotation of the solution (2.32) in the COV-gauge: 



where the gauge rotation U{x) is chosen such that A'^ — 0, i.e.. 



a = - [/t (d+U) 
9 



Eq. ( ^.37 ) is easily inverted to give 

[/^(a;~ , xj^) = P cxp < j dz^ a{z^ , x ±) 



(2.36) 



(2.37) 



(2.38) 



From eq. (2.3(;), is obtained indeed in the form (2.31), with U given in 



eq. (2.3S). The lower limit 



-cx) in the integral over x in eq. (2.38) 



has been chosen such as to impose the "retarded" boundary condition ( 2.20D . 
Furthermore: 



T+\x) = d+A'{x) = U{x)T+\x)U\x). 



(2.39) 



Together, eqs. (2.31), (2.34) and (2.38) provide an explicit expression for the 
LC-gauge solution A^ in terms of the colour source p in the COV-gauge. The 
corresponding expression in terms of the colour source in the LC-gauge p cannot 
be easily obtained: Eq. (2.33) implies indeed 



pU, 



(2.40) 



which implicitly determines a (and thus U) in terms of p, but which we don't 
know how to solve explicitly. But this is not a difficulty, as we argue now: 
Recall indeed that the classical source is just a "dumm y" v ariable which is 
integrated out in computing correlati ons according to eq. ( |2.6| ). Both the mea- 
sure and the weight function in eq. (2.6) are gauge invariant. Thus, one can 
compute correlation functions in the LC-gauge by performing a change of vari- 
ables p p, and thus replacing the a priori unknown functionals A^ [p] by the 
functionals -4' [/5], which are known explicitly. In other terms, one can replace 
eq. (P) by 



{A^^{x+,x)A''{x+,y)---)r = V[p] WAp]A^MA:[p] 



(2.41) 



where A'^lp] is the classical solution in some generic gauge (e.g., the LC-gauge), 
but expressed as a functional of the colour source p in the CO V-gange. 
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Moreover, the gauge-invariant observables can be expressed directly in terms of 
the gauge fields in the COV-gauge, although the corresponding expressions may 
look more complicated than in the LC-gauge. For instance, the operator which 
enters the gluon distribution can be written as (cf. eq. ( p.39| )) 

Tr {T+^{x)T+\y)} - Tr [u{x)T+\x)U\x) U{y)T+\y) UHy)] , (2.42) 

where the classical fields are in the LC-gauge in the l.h.s. and in the COV-gauge 



in the r.h.s, and U and W are given by eq. (2.38). Both writings express the 
gauge-invariant operator ( |2.18| ) (with the path 7 in Fig. ^ in the indicated 
gauges. (Indeed, U^{x,y) — W{x)U{y) for the COV-gauge field A'^ = 5^+a.) 
Note that, while in the LC-gauge the non-linear effects are encoded in the electric 
fields in the COV-gauge they are rather encoded in the Wilson lines U and 

(the corresponding field = —d^aa being linear in pa). 
Up to this point, the longitudinal structure of the source has been arbitrary: 
the solutions written above hold for any function p°(a;~). For what follows. 



however, it is useful to recall, from Sect. 2.1, that p has is localized near 
x^ = 0. More precisely, the quantum analysis in Sect. |3.4| will demonstrate 
that the classical source at the longitudinal scale has support at positive x~ , 



with < X < l/k~^- From eqs. ( 2.33 )-( 2.34 ), it is clear that this is also the 



longitudinal support of the "Coulomb field" a{x). Thus, integrals over x as 



that in eq. ( 2.38| ) receive contributions only from x in this limited range. The 



resulting longitudinal structure for the classical solution is illustrated in Fig. y_2|, 
and can be approximated as follows: 

A'ix-,x^) w 9{x-)-Vid'V^) = 0{x~)Al,{x^), (2.43) 

g 

T+\x) = d+A' « S{x~)Al,{x±). (2.44) 

It is here understood that the 6 and 9 functions of x~ are smeared over 

a distance Ax~ ~ l/fc"*". In the equations above, V and are the asymptotic 
values of the respective gauge rotations as x~ ^ 00 : 

VHx j_) = Feiip\ ig [ dz^ aiz^ ,xj_)\ . (2.45) 



In practice, U{x^,x±) — V{x±) for any x^ ^ l/fc+. Note that ( 2.45| ) is the 
same Wilson line as in the discussion of the eikonal approximation in Sect. 2.2.2 
(compare to eq. ( |2.25| ) there). In the present context, the eikonal approximation 



is implicit in the special geometry of the colour source in eq. (2.1), which is 
created by fast moving particles. 



2.4 The gluon distribution of the valence quarks 

To compute observables in the effective theory, one still needs an expression for 
the weight function Wr [p] ■ Before discussing the general construction of Wr [p] 
in Sect. 0, let us present a simple model for it, due to McLerran-Venugopalan 
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X 



Figure 12: The longitudinal structure of the colour source p and of the classical 
field solution for the effective theory at the scale A:"*" . As functions of cc" , a 
and are as localized as p. 



(MV) [P, which takes into account the colour charge of the valence quarks 
alone. That is, it ignores the quantum evolution of the colour sources with r. 
This model is expected to work better for a large nucleus, with atomic number 
A ^ 1 ; indeed, this has many valence quarks {A x Nc), and therefore as many 
colour sources, which can create a strong colour field already at moderate values 
of X, where the quantum evolution can be still neglected. In this model, t is 
fixed, but one can study the strong field effects (in particular, gluon saturation) 
in the limit where A is large. Besides, the MV model provides a reasonable 
initial condition for the quantum evolution towards small x, to be described 
later. 

The main assumption of the MV model is that the A x Nc valence quarks can 
be treated as independent colour sources. This relies on confinement. Note first 
that confinement plays no role for the dynamics in the transverse plane: Indeed, 
we probe the nucleus with large transerse momenta ^qcd^ ^^^^ is> O'^^'^ 

distance scales much shorter than those where confinement sets in. On the other 
hand, even at moderate values of x, we are still probing an integrated version 
of the hadron in the longitudinal direction, i.e., we measure all the "partons" 
(here, valence quarks) in a tube of transverse area AS± ~ 1 /Q^ and longitudinal 
extent Aa;^ ~ l/xP+ > 1/P+. The number of valence quarks which are crossed 
by this tube, 

AN^n AS± = AS± ^ ~ A^^^ , (2.46) 
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(with n = the number of quarks per unit transverse area, R the radius of a single 
nucleon, and Ra = A^/'^R the radius of the nucleus) increases with A, but these 
quarks are confined within different nucleons, so they are uncorrelated. When 
the number of partons A7V is large enough, the external probe "sees" them as a 
classical colour source with a random distribution over the transverse area. The 
total colour charge Q° in the tube is the incoherent sum of the colour charges 
of the individual partons. Thus, 



= 0, = g^CfAN = ASi_ (2.47) 



where we have used the fact that the colour charge squared of a single quark is 
g^^a^a _ g'^Cf. One can treat this charge as classical since, when AA^ is large 
enough, we can ignore commutators of charges: 

I [Q«, Q''] I = I ^r''^Q= I « Q\ (2.48) 

In order to take the continuum limit (i.e., the limit where the transverse area 
AS± of the tube is small^), it is convenient to introduce the colour charge 



densities p"(a; ,x±) (with the same meaning as in Sect. 2.1) and 

p'^ix^)- J dx-p°'{x-,x±) (2.49) 

(the colour charge per unit area in the transverse plane). Then, 

Q'* / d^x^p''{x±)= [ d^x± [ dx- p''{x-,x±), (2.50) 
Jas± Jas± J 



and eqs. ( ^.47D imply (recall that C/ = {N^ - l)/2Nc) 



g^A 

{Pa{xA_)pb{yA_))A = SabS''^\xj_ - y±_) PA, PA = o- n2 ' 



27ri?^ 



{pa{x ,x±)pb{y ,y±))A = SabS'^^''{x± - y±)S{x -y ) XAix ), 

dx^ Xa{x^) = PA- (2-51) 



Here, pA ~ A^^^ is the average colour charge squared of the valence quarks per 
unit transverse area and per colour, and Xa{x~) is the corresponding density 
per unit volume. The latter has some dependence upon x^ , whose precise 
form is, however, not important since the final formulae will involve only the 
integrated density pA- There is no explicit dependence upon x± in pA or Xa{x~) 
since we assume transverse homogeneity within the nuclear disk of radius Ra- 

^This amounts to increasing Q^, so, strictly speaking, at this step one should also include 
the DGLAP quantum evolution (i.e., the fact that, with increasing transverse resolution, the 
original "quark" is resolved into a set of smaller constituents). The quantum an alysis to be 



discussed later will include that in the "double-log approximation"; see Sect. 5.3 
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Finally, the correlations are local in x~ since, as argued before, colour sources 
at different values of x~ belong to different nucleons, so they are uncorrelated. 
All the higher-point, connected, correlation functions of Pa{x) are assumed to 
vanish. The non-zero correlators ( p. 51 ) are generated by the following weight 



function 



which is a Gaussian in pa, with a local kernel. This is gauge-invariant, so the 
variable pa in this expression can be the colour source in any gauge. The integral 



over xj_ in eq. (2.52) is effectively cutoff at Ra- By using this weight function, 
we shall now compute the observables introduced in Sect. ^.2[ . 
Consider first the gluon distribution in the low density regime, i.e., when the 
atomic number A is not too high, so that the corresponding classical field is 
weak and can be computed in the linear approximation. By expanding the 



general solution (2.31) to linear order in p, or, equivalently, by directly solving 



the linearized version of eq. ( p.l| ), one easily obtains 

k++ie kl ' ^ ' ~ ^kl 



which together with eq. ( 2.51 ) implies 



{K+{k)Tl+{~k))A ^ ^ {paik)Pa{-k))A = 7rR\{N! - 1)^ . (2.54) 



By inserting this approximation in eqs. ( 2.23| ) and ( |2.16| ), one obtains the fol- 



lowing estimates for the gluon density and distribution function: 

AfAik^) - (2.55) 
xG(x,Q^) . M^,,/"^ dl^^a^AN^^^ 



4^ Ja%^„ kl A^cD ' 

(with as — (7^/47r). The integral over k± in the second line has a logarithmic 
infrared divergence which has been cut by hand at the scale Aqcd since we 
know that, because of confinement, there cannot be gluon modes with transverse 
wavelengths larger t han 1 /Aqcd (see also Ref. Q). 

We recognize in eq. ( 2.55| ) the standard bremsstrahlung spectrum of soft "pho- 



tons" radiated by fast moving charges In deriving this result, we have 
however neglected the non-Abelian nature of the radiated fields, i.e., the fact 
that they represent gluons, and not photons. This will be corrected in the next 
subsection. 



2.5 Gluon saturation in a large nucleus 

According to eq. ( ^.55 ), the gluon density in the transverse phase-space is pro- 
portional to A^/"^, and becomes arbitrarily large when A increases. This is 



28 



however an artifact of our previous approximations which have neglected the 
interactions among the radiated gluons, i.e., the non-linear effects in the classi- 
cal field equations. To see this, one needs to recompute the gluon distribution 
by using the exact, non-linear solution for the classical field, as obtained in Sect. 



2.3, This involves the following LC-gauge field- field correlator 



yl A 



(2.56) 



which, in view of the non-linear calculation, has been rewritten in terms of 
the classi cal fi eld in the COV-gauge (cf. eq. ( ^.42 )), where T^"^ = —d^Ua- To 
evaluate (2.56), one expands the Wilson lines in powers of a and then contracts 
the a fields in all the possible ways with the following propagator: 



{aa{x)ab{y)) A = SabSix -y )-fAix ,x±-y±), 



7a(x ,k±) = ^a{x ). 



(2.57) 



We have used here p°'{x ,k±) — k'j_a°'{x ,k±), cf. eq. (2.34), together with 
eq. ( |2^ ) which holds in any gauge and, in particular, in the COV-gauge. The 
propagator ( ^.57 ) is very singular as fcj^ 0, but this turns out to be (almost) 
harmless for the considerations to follow. 

The fact that the fields a are uncorrelated in x~ greatly simplifies the calculation 
of the correlator ( 2.561 ). Indeed, this implies that the two COV-gauge electric 
fields d''ab{x) and d^ady) can be contracted only together, and not with the 
other fields a generated when expanding the Wilson lines. That is: 



(t^aV'"') _ (t^ac^'" j^.) = {d'^\x)d^c,dy)) lul{x)UUy) 
= 5{x- - y-)(Tr U\x)U{y)) ( - VI^a{x- ,x^ - y^)) , 



(2. 



where we have used — Uca in the adjoint representation. Eq. (2.58) can be 



proven as follows: i) By rotational symmetry, d^a{x) cannot be contracted with 
a field a{z~,x±) resulting from the expansion of U\x); indeed: 



{a{z ,x±)d^a{x ,x±)) cc d^jAix ,r±) 



oc 



-ik' 



(27r)2 ki 



0. 



ii) Contractions of the type 

{a{z^ ,y±)d^a{x~ , x±)) {a{u^ ,x±)d''a{y^ ,y±)) oc S{x^ - z")(5(u" - y^) 

are not allowed by the ordering of the Wilson lines in x" : a{z~,yj_) has been 
generated by expanding C/^(y), which requires < y~ (and similarly <x^). 
Then, the first contraction in ( |2.59 ) implies x^ — z^ < y~ , while the second 
one leads to the contradictory requirement y~ = u~ < x^ . 
The allowed contractions in eq. (2.58) involve: 



Sa{x ,xj_-yj_) = j^^—j{TrU\x ,xj_)U{x ,y±))A, 



(2.59) 
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which is hke the S'-matrix element (2.26) for the dipole-hadron scattering, but 
now for a colour dipole in the adjoint representation (i.e., a dipole made of two 
gluons). This can be computed by expanding the Wilson lines, performing con- 



tractions with the help of eq. (2.57), and recognizing the result as the expansion 



of an ordinary exponential. One thus finds (see also Sect. 5.1 for a more rapid 
derivation) : 



exp 



dz- -iA{z',ri_), (2.60) 



Sa{x ,r_L) 

where the exponent can be easily understood: It arises as 
[gT^iaaix) - aa{y))gT^{ai,{x) - atiy)) 



(2.61) 



where igT°'(aa{x) — aa{y)) is the amplitude for the dipole scattering off the 
"Coulomb" field aa, to lowest order in this field (i.e., the amplitude for a single 
scattering). Then, ( p.6l[ ) is the amplitude times the complex conjugate am- 
plitude, that is, the cross section for such a single scattering. This appears 
as an exponent in eq. ( 2.60| ) since this equation resums multiple scatterings to 
all orders, and, in the eikonal approximation, the all-order result is simply the 
exponential of the lowest order result. Since, moreover, is the field created 
by the colour sources in the hadron (here, the valence quarks), we deduce that 
eq. ( 2.60| ) describes the multiple scattering of the colour dipole off these colour 
sources. 

If the field aa is slowly varying over the transverse size rj_ = xj_ — yj_ of the 
dipole ( "small dipole" ) , one can expand 

gT%aa{x) - aM) ~ -gT'^ix' - y'05*«a(x) = gT-{x' - y')Tt\x), (2.62) 



and then eq. (2.61) involves the correlator of two (COV-gauge) electric fields. 
This is indeed the case, at it can be seen by an analysis of the exponent in 
eq. (g) : 



iA{x ,Oj_)-^A(a; ,r_L) 
Ha{x^) 



flAix ) 



(Pkj_ 1 



(27r)2 kl 



(2.63) 



The above integral over k± is dominated by soft momenta, and has even a 
logarithmic divergence which reflects the lack of confinement in our model (see 
also [p5l). Note, however, that the dominant, quadratic, infrared divergence 
~ J {(rk±/k']_), which would characterize the scattering of a co/ourerf particle (a 
single gluon) off the hadronic field^, has cancelled between the two components 

^Such a divergence would occur in {U^{x))a, which describes the scattering of a single 
gluon. 
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of the colourless dipole. The remaining, logarithmic, divergence can be cut off 
by hand, by introducing an infrared cutoff Aqcd- Then one can expand: 



(Pkj_ 1 - e 



ik\ -ri 



l/rl 



(27r)2 



fci 



(27r)2 ■fcj 2 



Idtt 



,2 A2 



(2.64) 



(This is valid to leading logarithmic accuracy, since the terms neglected in this 
way are not enhanced by a large transverse logarithm.) We thus obtain: 



SAix ,rj_) 



exp ■ 



— - — rj^^iA(x ) In- 



(2.65) 



which together with eq. (|2.58|) can be used to finally evaluate the gluon density 
( 2.23 ). This requires a double Fourier transform (to fc+ and k±), as shown 
in eq. ( 2.17 ). The presence of the (5- function in eq. ( 2.58 ) makes the Fourier 
transform to fc+ trivial, and one gets: 



J^Aik±) = 



47^3 



(fr±e 



dx SAix ,rj_){~\/j_-/A{x ,rj_)), 



(2.66) 



where (cf. eqs. ( ^.57 ) and ( 2.63 )) 



^1ia{x ,rj_) = Xa{x ) 



d'^p± e*' 
(2^~ 



Pi 



= — In ■ 



4^ rl^QCD 



dp.A{x ) 



(2.67) 



The non-linear effects in eq. are encoded in the quantity S a{x 

which finds its origin in the gauge rotati ons i n the r.h.s. of eq. ( 2.56|) . In 
fact, by replacing Sa{x~ jr±_) — > 1 in eq. (2.66), one would recover the linear 
approximation of eq. (2.55). To perform the integral over x" in eq. (^!66| ), 
we note that the quantity (2.67) is essentially the derivative w.r.t. x~ of the 
exponent in Sa{x^ ^r^^, eq. (2.65). Therefore: 



^fA{kl_] 

where 



Ar2_i /• l-exp| - ir^^g^ln^^-^^^ | 

47r4 i''"^' a,iV,ri 



asNcfiA = a^Nc / dx Xa{x ) - A^/^. 



(2.68) 



(2.69) 



Eq. (2.68) is the complet result for the gluon density of a large nucleus in the 
MV model To study its dependence upon fcj^, one must still perform 

the Fourier transform, but the result can be easily anticipated: 
i) At high momenta k± 3> Qa, the integral is dominated by small distances 
r±<til/QA, and can be evaluated by expanding out the exponential. To lowest 
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non-trivial order (which corresponds to the hnear approximation), one obtains 
the bremsstrahlung spectrum of eq. (2.55): 



UA{k±) oc 



kl 



for fc± > Qa- 



(2.70) 



ii) At small momenta, k± <^Qa, the dominant contribution comes from large 
distances r± ^ 1/Qa, where one can simply neglect the exponential in the 
numerator and recognize 1 /r^ as the Fourier transform^ of In fc^ : 



^fA{k±) 



— In for fc_L < Qa- 

asNc k\^ 



47r3 



(2.71) 



There are two fundamental differences between eqs. ( |2J0| ) and ( |2Jll ), which 
refer both to a saturation of the increase of the gluon density: either with 1/k\ 
(at fixed atomic number A)^ or with A (at fixed transverse momentum fcj^). In 
both cases, this saturation is only marginal : in the low-fc_L regime, eq. ( |2.7lD , 
the gluon density keeps increasing with and also with A, but this increase 

is only logarithmic^ in contrast to the strong, power-like, increase cx [A^/'^ /kY) 
in the high-A:j_ regime, eq. (2.70). 

Moreover, the gluon density at low k±^ is of order l/a^ , which is the maximum 
density allowed by the repulsive interactions between the strong colour fields 
— ^ lyA^ A^) ^ 1/ g. When increasing the atomic number A, the new gluons 
are produced preponderently at large transverse momenta ^ Qa, where this 
repulsion is less important. This is illustrated in Fig. |l^. 

To be more precise, the true scale which separates between the two regimes 
( p. 70 ) and ( p. 71 ) is not Qa, but rather the saturation momenturnQJA) which 
is the reciprocal of the distance l/r±_ where the exponent in eq. ( 2.68| ) becomes 
of order one. Thus, this is defined as the solution to the following equation: 



Ql{A) = - asN,^,A ln§^ 



(2.72) 



To clarify its physical interpretation, note that, at short-distances r±<^l/QA, 



fiA In ■ 



1 xG(x,l/r2^) 



(2.73) 



is the number of gluons (of each colour) having tranverse size r± per unit of 
transverse area (cf. eq. (p. 55 )). Since each such a gluon carries a colour charge 



squared {gT''){gT" 



g Nc, we deduce that 
o^sNct^A In- 



,2 A 2 



^The saturation scale provides tlie ultraviolet cutoff for the logari 



the short distances r± <^ 1/Qa a-r^ cut off by the exponential in eq. (2.68). 



(2.74) 



hm in cq. ([2.7l[) since 
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Figure 13: The gluon phase-space density AfA{k±) of a large nucleus (as de- 
scribed by the MV model) plotted as a function of k± . 



is the average colour charge squared of the gluons having tranverse size r± per 
unit area and per colour. Then, eq. (2.72) is the condition that the total colour 
charge squared within the area occupied by each gluon is of order one. This is 
the original criterion of saturation by Gribov, Levin and Ryskin , for which 
the MV model offers an explicit realization. 

To conclude this discussion of the MV model, note that, in the previous compu- 
tation, we have also obtained the ^'-matrix element SA{r±) for the dipole-hadron 
scattering (cf. Sect. p.2.2| ). This is given by eq. ( 2.65| ) with ^^(a;^) fiA and 
Nc = T'^T'^ replaced in general by the colour Casimir t"^t°- for the representation 
of interest (e.g., Cj = (N^ — l)/2Nc for the fundamental representation). As 



discussed after eq. ( 2.61 ), this describes the multiple scattering of the colour 
dipole on the col our fi eld in the hadron (here, the field of the valence quarks). 
According to eq. ( 2.65| ), one can distinguish, here too, between a short-distance 
and a large-distance regime, which moreover are separated by the same "satu- 
ration scale" as for the gluon distribution: 

i) A small-size dipole r±^l/Qs is only weakly interacting with the hadron: 



,2 \ 2 



for r_L < l/Qs{A), 



(2.75) 



a phenomenon usually referred to as "colour transparency" . 

ii) A relatively large dipole, with r±^l/Qs, is strongly absorbed: 



SAir±) 







for rj_ > l/QsiA), 



(2.76) 
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a situation commonly referred to as the "black disk", or "unitarity", limit. 
The remarkable fact that the critical dipole size is set by the saturation scale 
Qs can be understood as follows: A small dipole — small as compared to the 
typical variation scale of the external Coulomb field — couples to the associated 
electric field (cf. eq. (2.62)), so its cross-section for one scattering, eq. (2.61), 
is proportional to the number of gluons {J^~^''J-~^^) within the transverse area 
r'^ explored by the dipole. This is manifest on eq. (2.65), whose exponent is 
precisely the c olour charge squared of the gluons within that area (cf. the remark 
after eq. ( 2.74 )). At saturation, this charge becomes of order one, and the dipole 
is strongly interacting. The important lesson is that the unitarity limit (2.76) 
for the scattering of a small dipole on a high energy hadron is equivalent to 
gluon saturation in the hadron wavefunction |^, ^, Q . 



3 Quantum evolution and 

the Colour Glass Condensate 

In this section, we show that the classical Yang-Mills theory described in Sect. 
H can be actually derived from QCD as an effective theory at small x. This 
requires integrating out quantum fluctuations in layers of p^, which can be 
done with the help of a renormalization group equation (RGE) for the weight 
function VFt[p]- We shall not present all the calculations leading to this RGE; 
this would require heavy technical developments going far beyond the purpose 
of these lectures. (See Ref. Q for more details.) Rather, we shall emphasize 
the general strategy of this construction and the physical picture behind it (that 
of the colour glass), together with those elements of the calculation which are 
important to understand the structure of the final equation. 



3.1 The BFKL cascade 



In Sect. 2.1, we have argued that the radiation of a soft gluon by a fast parton 
via the tree-level graph shown in Fig. ^.a can be described as a classical 
process with a colour source whose structure is largerly fixed by the kinematics. 
Our main goal in this section will be to show that this picture is not spoilt 
by quantum corrections. Wc start by showing that the dominant quantum 
corrections, those which will be resummed in what follows, preserve indeed the 
separation of scales which lies at the basis of the effective theory developed in 
Sect. |. 

Consider first the lowest-order radiative correction to the tree-level graph in 
Fig. |lO|.a, namely, the emission of one additional (quantum) gluon, as shown 
in Fig. [l4|.a. At the same level of accuracy, one should include also the vertex 
and self-energy corrections illustrated i n Fi g. [T^.b, c. This will be done in the 



complete calculation presented in Sect. 3.4. But in order to get a simple order- 



of-magnitude estimate for the quantum corrections — which is our purpose in 
this subsection — it is enough to consider the radiative process in Fig. |l4[a. 
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Figure 14: Lowest-order quantum corrections to the emission of a soft gluon by 
a fast parton: a) a real-gluon emission; b) a vertex correction; c) a self-energy 
correction. 



The probabiUty for the emission of a quantum gluon with longitudinal momen- 
tum in the range > Pi > is 

AP cx / = In - a,ln-. 3.1 

TT p+ TT k+ X 

This becomes large when the available interval of rapidity Ar — ln(l/x) is large. 
This is the typical kind of quantum correction that we would like to resum here. 
A calculation which includes effects of order {as ln(l/x))" to all orders in n is 

said to be valid to "leading logarithmic accuracy" (LLA). 

The typical contributions to the logarithmic integration in eq. ( |3.l[ ) come from 
modes with momenta p^ deeply inside the strip: p"*" 3> p^ ^ ■ Thus, in 



Fig. 14, a, the soft final gluon with momentum fc+ is emitted typically from a 
relatively fast gluon, with momentum p^ ^ fc+. This latter gluon can therefore 
be seen as a component of the effective colour source at the soft scale k~^. In 
other terms, one can visualise the combined effect of the tree-level process, Fig. 
|lO|.a, and the first-order radiative correction. Fig. |l^.a, as the generation of a 
modified colour source at the scale A;"'', which receives contributions only from 
the modes with longitudinal momenta much larger than fc+. This is illustrated 
in Fig. |l^. 





Figure 15: Effective colour source after including the lowest-order radiative 
correction. 



Clearly, when x is small enough, ln(l/x) ~ l/ctg, the "correction" (3.1) becomes 
of 0(1), and it is highly probable that more gluons will be emitted along the 
way. This gives birth to the gluon cascade depicted in Fig. |l^.b, whose domi- 
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nant contribution, for a fixed number of "rungs" N, is of order (as ln(l/x))^, 
and comes from the kinematical domain where the longitudinal momenta are 
strongly ordered: 

p+ = p+ > p+ > p+ > • • • > p+ = fc+. (3.2) 

(Other momentum orderings give contributions which are suppressed by, at 
least, one factor of l/ln(l/x), and thus can be neglected to LLA.) With this 
ordering, this is the famous BFKL cascade, that we would like to include in our 



effective source. This should be possible since the hierarchy of scales in eq. (3.2) 
is indeed consistent with the kinemati cal a ssumptions in Sect. |[ 



Note first that, the strong ordering (3.2) in longitudinal momenta implies a 



corresponding ordering in the lifetimes of the emitted gluons (cf. eq. (2.2)): 

Ax^ > Axi > Ax+ > • • • > Ax^. (3.3) 

Because of this, any newly emitted gluon lives too shortly to notice the dynamics 
of the gluons above it. This is true in particular for the last emitted gluon, with 
momentum k~^, which "sees" the N previous gluons in the cascade as a frozen 
colour charge distribution, with an average colour charge Q = \J (QaQa) ~ 
Thus, this (A^+l)th gluon is emitted coherently off the colour charge fluctuations 
of the N previous ones, with a differential probability (compare to eq. ( ^.l| )) : 

dPN oc ^^N{T)dT. (3.4) 

TT 

When increasing the rapidity by one more step, t t + dr, the number of 
radiated gluons changes according to 

N{t + dr) = (1 + N{T))dPN + N(t){1 - dPN), (3.5) 

which together with eq. (3.4) implies (with = asNc/i:) 



xG(x,Q^) = ~ Case""""^. (3.6) 
dr 

Thus, the gluon distribution grows exponentially with t = ln(l/x). A more 
refined treatment, using the BFKL equation, gives k — 4 In 2, and shows that 



the prefactor C in the r.h.s. of eq. (3.6) has actually a weak dependence on r : 
C(x(a.r)-i/2 §1. 

Thus, the BFKL picture is that of an unstable growth of the colour charge 
fluctuations as x becomes smaller and smaller. However, this evolution assumes 
the radiated gluons to behave as free particles, so it ceases to be valid at very 
low X, where the gluon density becomes so large that their mutual interactions 
cannot be neglected anylonger. This happens, typically, when the interaction 



probability for the radiated gluons becomes of order one, cf. eq. (1.12), which 
is also the criterion for the saturation effects to be important (compare in this 
respect eq. ( 1.12| ) and eqs. (2.72)-( p.73| )). Thus one cannot study saturation 



consistently without including non-linear effects in the quantum evolution. It is 
our main objective in what follows to explain how to do that. 
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3.2 The quantum effective theory 

To the accuracy of interest, quantum corrections can be incorporated in the 
effective theory by renormahzing the source Pa and its correlation functions 
(i.e., the weight function Wt-[p]). The argument proceeds by induction: We 
assume the effective theory to exist at some scale A+ and show that it can be 
extended at the lower scale 6A+ ^ A+. Specifically: 

I) We assume that a quantum effective theory exists at some original scale A+ 
with A"*" <C ■ That is, we assume that the fast quantum modes with momenta 
p+ ^ A+ can be replaced, as far as their effects on the correlation functions 
at the scale A"*" are concerned, by a classical random source Pa with weight 
function Wa+ [p] . (We shall eventually convert A+ into the rapidity t by using 
T = ln(P+/A+).) On the other hand, the soft gluons, with momenta p+ < A+, 
are still explicitely present in the theory, as quantum gauge fields. Thus, this 
effective theory includes both the classical field A^[p\ generated by p, and the 
soft quantum gluons. 

Within this theory, the correlation functions of the soft (fc+ < A+) fields are 
obtained as (e.g., for the 2-point fimction) 



(TA^{x)A%y)) ^ VpWA+ip] 



J^* VA S{A+ ) At" {x)A'' (y) e *^[^^ "1 
J^^ VAS{A+) e'SlA,p] 



(3.7) 



where T stays for time ordering (i.e. ordering in a;+). This is written in the 

LC-gauge = 0, and involves two functional integrals: 

a) a quantum path integral over the soft gluon fields A^ at fixed p: 



, i VA5{A+)A^'{x)A''{y)e'S^'^^P^ , , 

{TA>^{x)A'^{y))p = ^ ^ ' ^ , 3.8 

/ VA5{A+) e^s[A,p] 



+ 



b) a classical average over p, like in eq. (2.6) 



{TA^^{x)A'^{y)) = VpW^^[p]{i:A^{x)A'^{y))p. (3.9) 



The upper script "A+" on the quantum path integral is to recall the restriction 
to soft (Ip^I < A+) longitudinal moment ap|. The action S[A, p\ is chosen such as 



to generate the classical field equations (2J_) in the saddle point approximation 



SS/SA^^ = 0. This requirement, together with gauge symmetry and the eikonal 
approximation, single out the following action ||3^ : 

S[A,p] = -ld^^\ F;,^Fr + d'xTT {p{x) W[A-]{x)}, (3.10) 

*The separation between fast and soft degrees of freedom according to their longitudi- 
nal momenta has a gauge-invariant meaning (within the LC-gauge) since the residual gauge 
transformations, being independent of x~ , cannot change the p+ momenta. 
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where W[A ] is a Wilson line in the temporal direction: 



W[A-]{x) ^Texpl^igJ da:+A~(a;)|. (3.11) 



With this action, the condition SS/SA^^ = implies indeed eq. ( |2.l| ) for field 
configur atio ns having A~ = 0. Thus, the classical solution A'^ = S^'^Al^[p] found 
in Sect. ^.3| is the tree-level field in the present quantum theory. 
As long as we are interested in correlation functions at the scale A"*", or slightly 
below it, we can satisfy ourselves with this classical (or saddle point) approxima- 



tion. That is, to the accuracy to which holds the effective theory in eq. (3.7), the 
gluon correlations at the scale A"*" can be computed from the classical field solu- 
tion, as in eq. (|2.6|). But quantum corrections become important when we con- 



sider correlations at a much softer scale ^ A"*", such that ag ln(A"'"/fc+) ~ 1. 

II) Within the quantum effective theory, we integrate out the semi-fast quantum 
fluctuations, i.e., the fields with longitudinal momenta inside the strip: 

6A+ < |p+| < A+, with 5<1 and a,. ln(l/fe) < 1. (3.12) 

This generates quantum corrections to the correlation functions at the softer 
scale bA~^, which can be computed by decomposing the total gluon field as 
follows: 

A>^ ^ A!i[p]+a>^ + SAi^. (3.13) 

Here, A'^ is the tree-level field, are the semi-fast fluctuations to be integrated 
out, and are the soft modes with momenta < 5A+ whose correlations 
receive quantum corrections from the semi-fast gluons. 

These induced correlations must be computed to leading order in asln(l/&) 
(LLA) , but to all orders in the classical fields A^ [p] (since we expect A^ ^ 1/g 
at saturation). This amounts to an one- loop calculation, but with the exact 
background field propagator {a^{x)a'^ [y]) p of the semi-fast gluons. For instance, 
the quantum corrections to the 2-point function read schematically: 

{{A^[p\ + 5A%A=[p\ + 5A^))^ - A^[p\A^[p] - 

= A'[p\{5A^)p + {5A')pA^[p\ + {5A'6A^)p (3.14) 

where the brackets (■ ■ •)p stand for the quantum average over the semi- fast fields 



in the background of p; this average is defined as in eq. (3.8), but with the 
functional integral now restricted to the fields a^. The purpose of the quantum 
calculation is to provide explicit expressions for the 1-point function {5A^) p and 
the 2-point function {5A^5A')p as functionals of p (to the indicated accuracy). 
Once these expressions are known, the 2-point function {A''{x)A^ {y)) at the 
scale &A+ can be finally computed as: 

{A^A^) = l^{{A}[p\ + SA^MA^P] + 5A^))p )^^, (3.15) 
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where the external brackets {■ ■ ■)wa denote the classical average over p with 
weight function Wa[p], as in eq. (|3.9[). 



Ill) We finally show that the induced correlations can be absorbed into a func- 
tional change W/^+ [p] W;,a+ [p] in the weight function for p. That is, the result 
( p. 15 ) of the classical+quantum calculation in the effective theory at the scale 



can be reproduced by a purely classical calculation, but with a modified 
weight function WbA+[p], corresponding to a new effective theory: 

(({A'[p] + SA^)iAnp] + = {A[pmp])^^^, (3.16) 



where the average in the r.h.s. is defined as in eq. ( |2.6| ), or (3.9), but with weight 
function T4^f,A+ [p] . This demonstrates the existence of the effective theory at the 
softer scale 6A+. 

Since AW = Wf,A+ — W\+ oc asln(l/6), the evolution of the weight function 
is best written in terms of rapidity: WV+Ar ^ W^r = —ArHWr, where r = 
ln(P"'"/A+), At = ln(l/6), and H = H[p, j^] is a functional differential operator 
acting on WV (generally, a non- linear functional of p). In the limit At 0, this 
gives a renormalization group equation (RGE) describing the flow of the weight 
function with r I03|, |36|l : 



dWrjp] 

dT 



-H 



5 

Sp 



Wr[p]. (3.17) 



By integrating this equation with initial conditions at r ^ 1 (i.e., at A"*" ~ 
P"*"), one can obtain the weight function at the rapidity r of interest. The 
initial conditions are not really perturbative, but one can rely on some non- 
perturbative model, like the MV model discussed in Sects. ^.4[-p^. 
A key ingredient in this approach, which makes the difference w.r.t. the BFKL 
equation, are the non-linear effects encoded in the background field calculation. 
Recall that p, and therefore the classical fields are random variables whose 

correlators ( |2.6| ) reproduce the gluon density and, more generally, the n-point 
correlation functions of the gluon fields at the scale A+. Thus by computing 
quantum corrections in the presence of these background fields, and then aver- 
aging over the latter, one is effectively studying quantum evolution in a medium 
with high gluon density. After each step in this evolution, the properties of the 
medium (i.e., the correlators of p) are updated, by including the latest quantum 
corrections. In terms of Feynman graphs of the ordinary perturbation theory, 
this corresponds to a complicated resummation of diagrams describing the inter- 
actions between the gluons radiated in different parton cascades and at different 
rapidities. A typical such a diagram is shown in Fig. h6l At low density, where 



the non-linear effects can be neglecte d, eg . ( 3.17 ) correctly reproduces the BFKL 



equation |3q], as it should (see Sect. |3.5| below). 
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Figure 16: A typical Feynman diagram that is imphcitely resummed in the 
quantum evolution of the effective theory. 



3.3 The Colour Glass Condensate 

Note the special form of the average in eq. ( |3.7| ). This is not the same as 

JVp Wa[p] f''VAA^'{x)A^'{y) e^^^I^^^l 
JVpWa[p] J^VAe'SlA,p] 



(3.18) 



In eq. (3.18), both the colour source pa and the gauge fields Ai^ are dynamical 
variables that are summed over on the same footing. They are free to take on 
values which extremize the total "effective action" : 



Seff[A, p] = S[A, p] - iIuWaIp]. (3.19) 

By contrast, in eq. ( t3.7[) , the average over Af^ is taken at fixed p : the gauge 
fields can vary in response to p, but p cannot vary in response to the gauge 
fields. That is, p is not a dynamical variable, but rather an "external" source. 
Giving a colour charge distribution Pa{x) specifies a medium in which propagate 
the quantum gluons. But this medium is, by itself, random, so after performing 
the quantum analysis at fixed p, one must also perform an average over p. The 
reason for treating p and A^^ differently lies is the separation of scales in the 
problem: the changes in p happens on time scales much larger than the lifetime 
of the soft gluons. This situation is typical for amorphous materials called 
"glasses" . 

The prototype of such systems is a "spin glass" that is, a collection of 

magnetic impurities (the "spins" ) which are randomly distributed in a non- 
magnetic metal host. For instance, one can take the spins to sit on a regular 
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lattice with lattice sites . . ., and interaction Hamiltonian 

Hj[S] = - ^ J^JS,S,, (3.20) 

<»,i> 

(the sum runs over all pairs < i,j >, and the spins Si are allowed to take two 
values, +1 or —1), but let their interactions (the "link variables" Jij) to be 
random, with a Gaussian probability distribution, for simplicity: 



dP[J] = n d.hjP{.hj), P{J^, 



1 



<i,3> 



v/27rA, 



(3.21) 



Physically, this corresponds to the fact that the modifications in occur on 
time scales much larger than the time scales characterizing the dynamics of the 
spins (e.g., their thermalization when the system is brought in contact with 
a thermal bath). In practice, the J^'s are frozen into their fixed values by 
rapid cooling when the sample is prepared. This kind or rapid cooling is called 
"quenching" , and one says that the Jij 's are "quenched variables" , as opposed to 
the "dynamical variables", the spins Si. This procedure selects random values 



for the Jy 's, with the probability distribution (3.21). 

Thus, the spins thermalizc for a given set of "quenched variables", and for each 
such a set one can compute the thermal partition function and the free energy: 



Z[J] 



{S} 



F[J] = -TlnZ[J]. 



(3.22) 



But the Jij 's are themselves random, so the experimentally relevant quantity is 
the following average 

F ^ {F[J]) = j dP[J]F[J] = j dP[J] lnZ[J]. (3.23) 

Note that it is lnZ[J], not Z[J] itself, which should be averaged ("quenched 
average"). Similarly, (connected) correlation functions are generated by the 
free energy in the presence of a site-dependent external magnetic field: 



{S,S,) - {S,){S.j) = / dP[J] 



\nZ[J,h] 
dhidhj 



(3.24) 



with \ nZ\J , h] defined as in eq. (|3.22|) , but with Hj[S] Hj[S] - J2^ htSi. 
Eqs. (^3|)-(|3^) are the analogs of eq. ( |3.7| ) for the problem at hand: the 
colour source pa is our "quenched variable" , and the quantum average over the 
field s at fixed p, eq. (3.8), corresponds to the thermal average at fixed Ji^'s, 
eq. ( 3.22|). As in eq. ( |3.23 ), it is InZ, and not Z, which is effec tivel y averaged 
in eq. (3/7) (the average of Z would rather corresponds to eq. ( 3.18| )). In fact, 
the connected correlation functions of the soft gluons in the effective theory are 
obtained from the following generating functional: 



Vp Wiy[p] In 



VA5[A+) e 



2S[A,p]-2 l j-A 



(3.25) 
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which is the analog of eq. (3.23) with In Z[J] In Z[J, h] . (The external current 
in (3.25) is just a device to generate Green's functions via differentiations, 
and should not be confused with the physical source pa-) 

We are thus naturally led to interprete the small-x component of the hadron 
wavefunction as a glass, with the colour charge density playing the role of the 
spin for spin glasses. Thus, this is a colour glass. Unlike what happens for spin 
glasses, which may have a non-zero value for the average magnetization (Si) (at 
least locally, i.e., at a given site), the average colour charge must be zero. 



{pa{x)) = at any x, 



(3.26) 



by gauge symmetry. In practice, this is insured by the fact that we sum over 
all the possible configurations of Pa{x) with a gauge-invariant weight function. 
Let us however examine a particular configuration Pa{x) from this ensemble. 
We now argue that, at sufficiently small x (or large atomic number A), this 
configuration describes typically a Base condensate. 

This applies to the saturated modes, i.e., the modes with transverse momenta 
^QCD ^ fc_L <C Qs{t) and longitudinal momenta fc+ = xP+ ^ P+. As argued 
in Sect. 2.5, these modes are characterized by a high giuon number density in 
the transverse phase-space, AfT{k±) ~ I /as- (This prediction of the classical 
MV model remains valid after including the quantum evolution, as we shall see 
in Sect. 5.4 below.) Microscopically, these modes correspond to bosonic states 
with large occupation numbers ~ l/ag- Each such a state is a Bose condensate. 
More precisely, the general definition of a Bose condensate is that of a quantum 
state in which the Fock space annihilation operator a*(/c) (cf. eq. (2.13)), or, 
equivalently, the field operator A'^(x), takes on a non-zero expectation value. 
This situation may be characterized as the spontaneous generation of a classical 
field. Of course, this cannot happen for gluons in the vacuum, as it would violate 
gauge symmetry. And, in an absolute sense, this cannot happen i n a h adron 
neither, since the average colour charge vanishes there too (cf. eq. ( 3.26 )), and 
therefore so does the associated classical field: (-4^[p]) = 0. But in the hadron 
there are colour sources, and, as argued before, they can be even treated as 
a classical charge distribution which is frozen during the short lifetime of the 
small-x gluons. Thus, over such a short time scale (short as compared to the 
typical time scale for changes in the colour distribution), one effectively has a 
non-t rivia l class ical fi eld -4J,[/)]. At saturation, this field is typically strong (cf. 
eqs. (|2J^) and §M)) : 



1 1 



ttNc 

1 



cesr]^ 



\J a^r 



for r_L > , 



(3.27) 
(3.28) 



and its typical amplitude (3.2?) at large rj^ is even independent of the actual 
strength p ^ {PaPa) ^ Qs of the colour source. This can be thus characterized 
as a Bose condensate. 
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We thus see that it is the same fundamental separation in time scales which 
allows us to speak about both the colour glass and the Bose condensate, although 
these two concepts seem at a first sight contradictory: the notion of a "glass" 
makes explicit reference to the average over p, while the "condensate" rather 
refers to a specific realization of p, before averaging. 



3.4 The renormalization group equation 



As explained in Sect. 3.2, the quantum evolution of the effective theory is 



obtained by matching correlations computed in two ways: (a) via a classi- 
cal+quantum calculation in the effective theory at the scale A+, and (b) via 
a purely classical calculation within the effective theory at the scale 6A+. The 
quantum corrections that are included in this way are those generated by the 
coupling between the "semi- fast" gluons with p+ momenta in the strip ( [3.12| ) and 
the "soft" gluons with momenta < 6A+. To the accuracy of interest, 
it is sufficient to consider the eikonal coupling SA^Spa to the plus component 
SJ^ = Spa of the colour current of the semi-fast gluons. Indeed, these gluons 
are relatively fast moving in the direction, so 5pa is the large component of 
their current. 

The results of the matching can be summarized as follows: 
i) T o 0{ as ln(l/b)), the induced correlations of the transverse fields AJ, (see 
eq. ( 3.14 ) for an example) can be all related to the following 1-point and 2- 
point functions of dp (with At = ln(l/6)) : 



1 

a7 



(3.29) 



dx / dy (Spa{x^,x)Spb{x+,y))p, (3.30) 



where, as in eq. ( 3.14 ), (• • •)p denotes the average over semi-fast quantum fluc- 
tuations in the background of the tree-level source p. 

Thus the quantum evolution consists in adding new correlations a and x to p. 
ii) These new correlations can be included in the weight function Wr[jo] by 
allowing this to evolve with r according to the following RGE |pq, • 



dWr[p] _ 1 

dT " 2 5p^{xA^)Spb{yi_) 



[WrXt] 



5p%{x±_ 



(3.31) 



We use here compact notations where cr^ = cra{x±), Xxy = Xab{x±,y±), and 
repeated colour indices (and coordinates) are understood to be summed (inte- 
grated) over. The notation p'^{x±) will be explained later (see eq. ( |3.46 )). 
A complete proof of the statements above would require the lengthy analysis 
of Refs. |3^. But assuming them to be true, it is easy to understand the 



general structure of the RGE (3.31). Indeed, according to eqs. (3.29)-(B.30), the 
induced correlations that we need to take into account are (with the notations 
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of eq. (p5|)): 



U{Pa+ 5pa)xAPb + 5pb)y±) „) - (Pa{x±)pb{y±))^^^ 

D[p]Wr[p] At {aa{x±)pb{y±) + Pa{x±)<Tb{y±) + Xab{x±,y±)} 



(3.32) 



1 



Pa{xi_)pb{yi_), 



where the colour mdices c, d (the transverse coordinates u_l) in the last line 
are to be summed (integrated) over. After a few integrations by parts w.r.t. p, 
the last expression can be recast into the form: 



D 



[P] Pa{xi_)pb{yi_)AWr[p\, 



(3.33) 



with AWV[p] given by the finite-difference version of eq. (3.31). 

In eqs. (B.32)-(3.33), we have considered only correlators of tw;o-dimensional 

(or "integrated") charge densities, like 



Pa{xA_) 



dx Pa{x ,X±), 



(3.34) 



and similarly Spa{x±_). This is in agreement with eqs. (3.29)^(3.30), which show 
that only such integrated (over x~) quantum corrections are relevant to the or- 
der of interest, and is moreover physically intuitive: The soft gluons (fc"*" bA^) 
to which applies the effective theory are unable to discriminate the internal lon- 
gitudinal structure of their sources, which are localized in x^ over relatively 
short distances <C 1/&A+, because of their large momenta. Although es- 
sentially correct, this argument is a little too simplistic as shown by the fact 
that some of the quantities encountered before are in fact sensitive to the lon- 
gitudinal struct ure o f p (i.e., they are not simply functionals of the integrated 
charge density ( 3.34 )). A generic example i s the bac kgrou nd field .4'[/9], or any 
other quantity built with the Wilson lines ( 2.38 ) or ( 2.45 ). Such quantities are 
sensitive to the x~ dependence of p because of the path-ordering of the Wilson 
lines in x~ . The ordering is important since colour matrices p{x~) = pa{x~)T°' 
at different values of x~ do not commute with each other. This suggests that 
the correct way to think of an "integrated" version of the hadron (over x~) is 
in terms of Wilson lines — which take into account the colour precession in the 
colour field of the hadron, with the proper ordering of colour matrices — , and 
not of 2-dimensional charge densities like (3.34). This will be confirmed by the 
subsequent analysis of the quantum corrections. 



3.4.1 The quantum colour source 

For the purposes of the quantum calculation, it is useful to expand the action 
S[A, p] = S[A -I- a -I- SA, p] to quadratic order in the small fiuctuations a^, and 
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retain only their eikonal coupling to the component SA^ of the soft fields: 
S[A + 5A + a,p] « S[A + 5A,p] + la^ix)G-^{x,y)[p]a''iy)-6A-Spa, 



with 



G-J^\x,y)[p] 



S'S[A, p] 



SA'^{x)6A^,{y) 



(3.35) 
(3.36) 



Spaix) 



S^S 



6Aaix)6A-iy) 



2SAaix)6A'^,iy)6A^{z) 



(3.37) 



where it is understood that only the soft modes with fc+ foA+ are kept in the 
products of fields. 



The expansion (3.35) corresponds to a one-loop approximation for the soft cor- 
relation functions like {6A^)p and {6A^SA^)p (cf. eq. ( |3.14D ), but where the 
propagator iG'^'^{x,y) = (Ta''(a;)a'^(?/))p of the semi-fast gluons running along 
the loop is computed in the background of the tree- level field A^ [p] , by inverting 
the differential operator in eq. ( 3.36| ). 





a) 



Figure 17: Some typical Feynman diagrams for x s-nd a. The internal wavy 
lines are propagators of the semi- fast gluons; the external dotted lines carry soft 
momenta, and couple to the fields 5A~ . (a) A contribution to x- The external 
blobs denote insertions of the electric field ; the internal line with a blob 
denotes the background field propagator, (b) A contribution to a to linear order 
in p. The continuous line represents the source p. 

To gain some more intuition, we use as an example the contributions to 6pa 
coming from the Yang-Mills piece of the action, Sym = J d'^x{—F^^/4:) : 



6paix)\ 



2gr'^T+\x)al{x) + g {d+ al{x))al{x) . 



YM 



(3.38) 



The first term in the r.h.s., which is linear in a*, is the only one to contribute 



to X, eq. (3.3C), to leading order in ag. It generates the tree-like diagram in 
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Fig. |T^.a, where the internal Hne with a blob represents the background field 
propagator G^^{x,y) of the semi-fast gluons. Physically, Fig. ^a describes the 
emission of an on-shell (or "real") semi-fast gluon by the classical source. 
Since (a*) = 0, it is only the second, quadratic term in the r.h.s. of eq. (3.38) 
which contributes to a, eq. ( 3.291 ). In Fig. |l^.b we show such a contribution of 
lowest order in p . (This involves also vertices from the Wilson line piece of the 
action, eqs. ( 3.1C| )-(3.11).) Obviously, this represents a vertex correction to the 
tree- level emission in Fig. 10, a. 

The structures illustrated by Figs, ^.a and b are generic: x is the "real" correc- 
tion, whose iteration generates the gluon cascades; a is the "virtual" correction, 
which provides one-loop corrections to the emission vertices in these cascades. 
Both X a-iid a include terms non-linear in p which describe interactions among 
gluons at different rapidities in different cascades. In general, real and virtual 
corrections are related by gauge symmetry, and this is also the case for x 
a, as we shall discuss later. 

The diagrams contributing to a and x in the general case, together with their 
explicit evaluation, can be found in Ref. Here, we shall present only the 

final results of this calculation. 



3.4.2 The induced colour source and field 



For the reasons explained in Sect. 2.3, it is more convenient to work with the 
colour source pa in the covariant gauge. The corresponding weight function 



Wr[/5] obeys an evolution equation similar to (3.31), but with modified coeffi- 
cients a and whi ch ar e obtained from the LC-gauge coefficients a and x via 
the gauge rotation ( 2.36| ). In what follows, we shall give directly the final results 
for these COV-gauge quantities. 

Consider first the induced source bp"" = {6p°')p, that is, the correction to the 
average colour charge density generated by the polarization of the semi-fast 
gluons. After rotation to the COV-gauge, this reads |^ : 



Spaix) = Fa{x )(- 

where the "form factor" 



-Vil^a(xi)), 



(3.39) 



ibAx 



iAx 



Fa{x-) 



(3.40) 



specifies the longitudinal profile of Spa, while {V^ = V^x±), cf. eq. ( ^.45| )) 



1 



d^z^ 



(27r)2 (x^^z^y 



■Tt:(t°-V^V^ 



(3.41) 



contains the dependence upon the background field Ua (via the Wilson lines V 
and V^), t ogether with the transverse and colour structure oi Spa- By comparing 
eqs. ( 3.3£ ) and ( 2.33 ), we deduce that F{x~)v°-{x±) is the induced field in the 
COV-gauge, i.e., the quantum correction to the tree- level field aa- Since: 



dx F\ (x 



Ar, 



(3.42) 
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eqs. ( 3.29| ) and ( 3.39| ) immediately imply : 

aa{x^) = -Vll^aix^). (3.43) 

This is the coefficient of the virtual term in the RGE for WV [p] . 

But the longitudinal structure of Spa is also interesting. Eq. ( |3.39| ) shows that 

the induced source and field have typically support at^ 

1/A+ < X- < 1/(6A+). (3.44) 

Recall that Spa has been generated by integrating out quantum fluctuations in 
the strip 6A+ <C <C A"*". Thus, when integrating out quantum giuons in 
layers oi p'^, one builds the classical source p (or field a) in layers of , with 
a one-to-one correspondence between the x~ coordinate of a given layer and 
the p"*" momenta of the modes that have been integrated out to generate that 
layer. By induction, we deduce that Pa{x) (= the colour source generated by the 
quantum evolution down to A^) has support at < x~ ^ 1/A^, as anticipated 
in Sect. ^. This allows us to consider only positive values for x~ in what follows. 
To exploit this tight correspondence between and x~ , it is convenient to use 
the space-time rapidity y, 

y = ln(a;"/a;o ), x^ = l/P+, -oo < y < oo , (3.45) 

to indicate the longitudinal coordinate of a field. We shall set, e.g., 

Py{x±) = x^p'^{x^,x±) for x^=Xy=XQe^, 

dyp^ix^) = J dx- p''ix-,x^), (3.46) 

and similarly for the other fields (p, a, etc.). The previous discussion on the 
longitudinal structure can then be summarized as follows: 
The source Py{x±) generated by the quantum evolution from r' = up to r 
has support at y in the interval < y < r. When new quantum modes, with 
rapidities r' in the interval r < r' < t + At, are integrated out, the preexisting 
colour source at y < r is not changed, but some new contribution is added to 
it, in the rapidity bin t < y < t + At. Because of that, AW = Wt+At — Wr 
involves only the change in py within that last bin. In the continuum limit 
At 0, this generates the functional derivatives of WV with respect to Py at 



y = T, as shown in eq. ( 3.31 ). This clarifies the longitudinal structure of the 
RGE. 



Consider also the transverse and colour structure of the induced field (3.41). 
This can be understood by reference to Fig. |l^.b. The transverse kernel in 
eq. ( 3.4l| ) has been generated as: 



1 1 1 — 1 a;* — 

(2^ (xi_ - zi_Y ^ 2^{x^- z^Y 2^ (a^-L - z±? ' 

^Indeed, Fj\^{x~) « both for small x~ ^ 1/A+ (since in this case the two exponentials 
mutually cancel), and for large x~ ^ 1/6A+ (where the two exponentials are individually 
small). 
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where (compare with eqs. ( ^.34 ) and ( ^.53 )) 



1 



27r (a;j_ ~ z±y 



(fp± —ip^ 



Pi 



^ip±-{x±-z±) 



(3.48) 



is the propagator of the semi-fast gluon emitted by the sou rce p (recaU that 
!F~^^ « {ip> /p'^)p to hnear order). The two Wilson hnes in eq. ( 3.4l| ) account for 
the scattering of this semi-fast gluon off the background field at z± (this brings 
in a factor V^{z±) in the eikonal approximation), and for its gauge rotation by 
the class ical field A^{x) at x~ > 1/A+ (cf. eq. ( 3.44 )), which is a pure gauge 
(cf. eq. (HI)). 



3.4.3 The RGE in the a— representation 



Eq. (3.43) suggests that it may be technically simpler and physically more trans- 
parent to work directly with the classical field aa and the quantum corrections 
to it (like i^a), rather than with the colour source pa and the corresponding 
corrections (like (Tq). This point of view is also supported by the fact that the 
LC- gau ge fi eld and the related observables are primarily related to aa (cf. Sects . 
2.2 and |2.3|) , and reexpressing them in terms of p — with the help of eq. ( 2.34 ) 
— would introduce a dependence upon the unphysical infrared cutoff p. 
For these reasons, we prefer to work in the a~representation, in which observ- 
ables are expressed in terms of a, and the average is performed with the weight 
function WV[q;] = Wr[p = — V^a]. This satisfies the following RGE, which is 
obtained after a change of variables in eq. (3.31) : 



dWr[a] 
dr 



1 



2 ,Sa° (xi )5a\ (y^) ' " '"^^ Sa^ {x^_ ) 
where z/^ = i'"-{x±), cf. eq. ( 3.41 ), and rj^y = ri°-''{x±,y±), with 



,j''\x^,yj_) = J d'zidV(xi|-^|2±)r'(^±, 



Uj_) {Uj_\ 



(3.49) 



(3.50) 



It is thus sufficient to give the result for the "real correction" {SpaSpb)p directly 
in the a-representation (cf. eqs. ( 3.3C| ) and (3.50)). This reads p7[ |: 



tPz_i (x* — z*)(y* — z') 
nj (2^ {x^-z^)^y^-z^y 

x{l + V^Vy - V^V. ~ VlVyY ■ (3.51) 



The transverse and colour structure of 77 have the same pattern as discussed 
after eq. ( |3.53 ) in connection with v. 

The r.h.s. of eq. (3.49) involves functional derivatives w.r.t. the colour field 
a°(a;_L) at the end point y = r. When applied to the coefficients 77 and this 
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requires the corresponding derivatives of the Wilson lines V and V\ that we 
compute now. Note first that, since cty = for y > r, we can rewrite 

vHxj_) = Pe*f /"oo-iy^yC^i) = pe^s/Ioody^yC^i). (3.52) 
Therefore (with 6xij = 6^^'>{x± — yi_)): 

^ ^^^-^ ^ J^) = -^9S.yVix^)T . (3.53) 



A simple interpretation of the four terms in cq. ( 3.51 ) follows from the dual 
picture of the dipole-hadron scattering, in which the quantum evolution is put 
in the dipole wavefunction, and, more generally, in the Wilson line operators 
through which a generic external projectile scatters off the hadronic target p^ , 
H, |6[ |7[J|, H, 1^, (See also the lectures notes by Al Mueller in this 
volume |14|.) Recent analyses of the high energy scattering from this dual 
perspective have led to a set of coupled evolution equations for the correlation 
functions of Wilson lines, originally derived by Balitsky Q (see also |50t), 
and subsequently reformulated by Weigert pTf in a compact way, as a functional 
evolution equation for the generating functional of these correlation functions. 
It turns out that Weigert's equation is equivalent to the RGE ( 3.49D jsT 



which demonstrates the equivalence between the two descriptions — the target 
picture and the projectile picture — of the nonlinear evolution in QCD at small 



X. We shall say more on Balitsky's equations in Sect. 4.3 



3.5 Recovering the BFKL equation 

Before studying more general properties and consequences of the RGE in the 
next section, let us rapidly show that, in the weak field (or low density) limit, this 
equation reproduces the BFKL equation, as expected |Q. Eq. ( 3.32| ) implies 
the following evolution equation for the 2-point function {pp)r '■ 

-^{Pa{x±)pb{vl_))r^{(^a{xi_)pb{Vi_) + Pa{xi_)ui,{y±_) + Xab{x ^ , y ^_)) r (3.54) 

For a generic, strong, source p, the coefficients a and x &re non-linear in p to all 



orders, so the r.h.s. of eq. (3.54) involves n-point correlators (p(l)p(2) • • • p{n))r 
of arbitrarily high order^ n. But in the weak field limit, where a is linear in 
p and X is quadratic, this becomes a closed equation for the 2-point function, 
which coincides with the BFKL equation, as we show now. 
Specifically, consider the evolution equation for the following 2-point function: 

^r(fci) = {pa{kl_)pa{-kl_))r, (3.55) 

^''incidentally, this shows that the n-point functions of p do not form a convenient basis 
to study the non-linearities in the evolution. By contrast, the correlator s of the Wilson lines 
form a more convenient such a basis as we shall discuss in Sect. 4.S. 
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{pa{k±) is the Fourier transform of Pa{x±)), which according to eqs. (2.16) and 
( ^.54 ) represents the "unintegrated gluon distribution" : 



dxGjjc, kl) 
dlnkl 



(3.56) 



in the weak field regime. In this regime, one can also expand ~ 1 + iga{x±), 
and therefore 

Then, eq. (3.41) reduces to (with as = g'^/iir and k^aa{k±) = Pa{k±)) ■ 



For X one obtains similarly 

Xi°j(fc±,-fci) = 4asN,kl 



kl 



(2^)2 pUp^ - k^Y ■ 

d'^P± Pa{P±)Pa{-P±) 
(2^)2 pi(fc^-pi)2 



(3.57) 



(3.58) 



By inserting eqs. ( |3.57 ) and (3.58) into the evolution equation (3.54), and using 
( p. 55 ), one finally obtains: 



dpT{k\) 
dr 



kl 



Mr (pi) 



27r2 



d^P± 



kl 



piikj_ -p±y 



■ Prik^ 



(3.59) 



which coincides, as anticipated, with the BFKL equation ^, ^. The first term 
in the r.h.s., which here is generated by x*"""*: is the real BFKL kernel, while the 
second term, coming from a^^\ is the corresponding virtual kernel. 
Note finally that the BFKL approximation has been obtained by expanding the 
Wilson lines to linear order in ga"", thus, this is formally the same as the lowest 
order perturbative expansion of the RGE. 



4 A functional Fokker- Planck equation 



We now dispose of a powerful tool — the functional RGE ( |3.49|) — to construct 
the effective theory by integrating out quantum fluctuations in perturbation 
theory. Eq. ( [3.49 ) has a rich and elegant mathematical structure, to be described 
in Sects. 4.1 and 4.2. Then, in Sects. 4.3 and ||, we shall indicate two strategies 
to make use of this equation: 

i) One can use it to derive ordinary (i.e., non- functional) evolution equations for 
the correlation functions of interest, like we did for the 2-point function {pp)t 
in Sect. \i.5\ When specialized to correlation functions of the Wilson lines, this 
strategy leads to a system of equations originally derived by Balitsky Q . This 
will be discussed in Sect. 4.2. 
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A difficulty witfi tfiis approach is that it generally leads to coupled equations 
(the 2-point function is coupled to the 4-point one, etc.), so that one has to 
follow simultaneously the evolution of infinitely many correlators. Still, some 
progress has been done, by using functional techniques ||5^, and, especially, by 
recognizing that, in the large Nc limit, a closed equation can be written for the 
2-point function: this is the Kovchegov equation p5{ . 

ii) One can try and solve directly the functional RGE, with appropriate initial 
conditions. An exact but formal solution can be written in the form of a path 
integral This is well suited for lattice simulations in 2+1 dimensions. But 
approximate analytic solutions, which allow for a more direct physical insight, 
have been found as well |l^, |57|. These solutions will be described in Sect. ||. 



4.1 General properties and consequences of the RGE 



We start with a summary of the most important properties of the RGE ( 3.49D 



i) The coefficients rj and v are real quantities. Moreover, r/ is sym- 
metric: riab{x±,y±) — riba{y±, x±), and positive semi-definite. 

ii) The RGE preserves the normalization of the weight function: 

VaWr[a] = 1 at any r. (4.1) 



Ind eed, the r.h.s. of eq. (3.49) is a total derivative with respect to a. Thus, if 
eq. ([4.l|) is satisfied by the initial condition at tq , it remains true at any t > tq . 
Properties (i) and (ii) guarantee that the solution WT-[a] to the RGE has a 



meaningful probabilistic interpretation (cf. the discussion prior to eq. (2.4)). 
iii) The momentum rapidity r and the space-time rapidity y are iden- 
tified by the quantum evolution. That is, the field ay in the rapidity bin 
(y, y + dy) is generated by the quantum evolut ion from r = yuptor = y + dy. 
This follows from the discussion in Sect. 3.4.2| , and implies that the two rapidi- 



ties can be treated as only one variable, the "evolution time" . 
With this interpretation, the function {ay(a:j^) | — oo < y < oo} — which phys- 
ically represents the longitudinal profile of the 3-dimensional field a° (x^ , x± ) 



in units of rapidity (cf. eq. (B.46)) — is viewed as a trajectory in the functional 
space spanned by the 2-dimensional fields a°'{x±). Quantum evolution then 
appears as the progression of the "point" a°'{x±) along this trajectory. Thus, 
eq. ( |3.49| ) describes effectively a field theory in 2-1-1 dimensions (the transverse 
coordinates and the "evolution time"), which is however non-local in both x± 



and y (since the coefficients (3.41) and (3.51) of the RGE involve ay at all the 



"times" y < T, via the Wilson lines (3.52)) 



iv) The initial condition. Let the quantum evolution proceed from some 
original "time" tq up to the actual "time" r. The "trajectory" {ay{x±) | — oo < 
y < oo} can be decomposed into three pieces: a) The field ay at y < tq belongs 
to the initial conditions, b) The field ay at tq < y < t is generated by the 
quantum evolution, c) There is no field at all at larger y: ay = for any y > r. 
Thus: 

Wr[a] = S[a>]Wr[a<], (4.2) 



51 



where (tty ) is the function ay for y < r (respectively, y > r) : 

ay{x^) = 6{T-y)af{xA_) + 6'(y - r)a> (x^l), (4.3) 

and the (5-functional (5[q!>] should be understood with a discretization of the 
configuration space, as in eq. (12.5^: 



%>] ^ nnn'^K(^-L))- (4-4) 

y>r a x± 

Moreover, it can be shown |5^ that yVT-[a^] has the factorized structure: 

Wr[a<] - Wr^rAa\V^] WrM , (4.5) 
where Wt-q [oi] is the initial weight function at tq , and 

V^{xi_) EE Ve'sPLdyo^yix^) (4g) 



is the Wilson line built with the initial field. In eq. (4.5), it is understood that, 
in Wro) the field argument ay has support at y < tq, while in Wt, ro it has 
support at To < y < T. The "propagator" 'Wt,to from tq to t depends also upon 
the initial field at y < t q, but only in an integrated way, via the Wilson lines Vb 
and Vq . From eq. ( |4.5| ) we deduce that W,-^ [ajVo] — > 1 when t — > tq. 
The initial weight function cannot be obtained within the present formal- 
ism, but rather requires some model for the hadron wavcfunction at rapidity tq. 
It is convenient to choose a moderate value for tq — ln(l/xo), e.g., xq ~ 10^^. 
This xo is small enough for the LLA to apply, but still large enough for the non- 
linear effects to remain negligible. Then one can use initial conditions which 



are consistent with the standard, linear, evolution equations (cf. Sect. 5.2 be- 
low). Once a convenient value for xq has been chosen, one can always redefine 
T = ln(xo/x) so that the initial condition is formulated at tq — 0. With this 
choice, the field ay at positive rapidities y > is generated by the quantum 
evolution, while the field at negative rapidities y < must be specified by the 
initial condition. 



v) The Hamiltonian structure of the RGE. Eq. (3.49) can be rewritten 
as: 

= ^ IV +11 -u'^\Wr\ (4 7) 

A crucial property, with many consequences, is that the second term within 
the braces is actually zero. Indeed, the following relation holds between the 




coefficients of the RGE |47, 37 
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It is easy to prove this relation by using eq. (3.53) to act with 8/5a\{y±^ on 
riab{x±,y±), eq. ( |3.5lD . This yields, e.g.. 



cb 



cb 



(4.9) 



= 0, 



where both terms in the second line vanish because of the antisymmetry of the 
colour group generators in the adjoint representation (e.g., {T^)ab — 0). The 
only nonvanishing contribution is 

6 



{V^V^r = -^9S.y{T''V^V.) = tgS:.yTr(T'^V^V.) , (4.10) 



which reproduces indeed eq. ( [3.41 ) after integration over y±, since 
{x'-z'){y''-z') 1 



IC{xj_,yj_,ZA_) 



{x±-z±)^{y±-z±)^ {x±-z±_y 



for yj_^x±. (4.11) 



With eqs. (13) and dJ), the RGE can be brought into a Hamiltonian form 

dWrla] 



dr 

with the following Hamiltonian: 



<Py± 



Ixy 



iS 



(Pz± 
2tt 



(4.12) 



<Px±_ z* — 
27r (zj_ ~ x±y- 



(1 - V;^T4)a6 



(4.13) 



which is Hermitian (since rj'^y is real and symmetric) and positive semi-definite 
(since the "current" J*(zj^) is itself Hermitian). 

vi) The infrared and ultraviolet behaviours of the RGE. These are deter- 
mined by the kernel 77"'' (xj^, y±^ in the Hamiltoni an. In the infrared limit, where 
Zi_ is much larger than both x±^ and (see eq. ( 3.51 )), /C(xj^, yj^, zj^) ~ ^/zj 



and the ensuing integral {d z±/zj_) has a logarithmic infrared divergent 
Thus, there is potentially an IR problem in the RGE. This is not necessar- 
ily a real difhculty, since IR problems are expected to be absent only for the 
gauge-invariant observables. We shall see indeed, on specific examples, that 
the IR divergences cancel when the RGE is used to derive evolution equations 
for gauge-invariant quantities. This cancellation relies in a crucial way on the 
property (4.8). 

Coming now to the ultraviolet, or short-distance, behaviour, it is easy to see on 
eq. (3.51) that no UV problem is to be anticipated. For instance, the would-be 
linear pole of IC{xj_, y±, z±) at |z_l — x±\ — > is actually cancelled by the factor 
1 — V^Vx which vanishes in the same limit. 



^^This infrared behaviour is not modified by the z± dependence of the Wilson hnes since, 
e.g., (V^Vz) ^ as \zj_ — xj_ \ —* 00; cf. Sect. 5.2 below. 
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4.2 Quantum evolution as Brownian motion 

To clarify the probabilistic interpretation of the RGE ( |3.49 ), we start by recall- 



ing the simplest example of a stochastic process, namely the Brownian motion 
of a small particle in a viscous liquid and in the presence of some external force, 
like gravitation The particle is so small that it can feel the collisions with 
the molecules in the liquid; after each such a collision, the velocity of the particle 
changes randomly. And the liquid is so viscous that, after each collision, the 
particle enters immediately a constant velocity regime in which the friction force 
oc v'^ (with the velocity of the particle) is equilibrated by the random force 
due to collisions together with the external force F^{x). In these conditions, the 
particle executes a random walk whose description is necessary statistical. The 
relevant quantity is the probability density P{x, t) to find the particle at point 
X at time t. This is normalized as: 

(fxP{x,t) = 1, (4.14) 

and obeys an evolution equation of the diffusion type, known as the Fokker- 
Planck equation |^ : 

dP{x,t) _ o/^ .^ ^ 



dt dx'^dx'^ 



P{x,t) - —(^F\x)P{x,t)). (4.15) 



Here, D is the diffusion coefficient, which is a measure of the strength of the 
random force; for simplicity, we assume this to be a constant, i.e., independent 



of X or t. The solution to eq. (4.15) corresponding to some arbitrary initial 



condition P[x,to) can be written as 

P{x,t) = J d^xoP{x,t\xo,to)P{xo,to), (4.16) 

where P{x,t\xo,to) is the solution to ( |4.15| ) with the initial condition: 

P{x,to\xo,to) ^ d'^^\x-xo). (4.17) 

Physically, this is the probability density to find the particle at point x at time 
t knowing that it was at xq at time Iq. 

If — 0, this solution is immediately obtained by going to momentum space: 
The Fourier transform P{k,t) of P{x,to\xo,0) = P{x — XQ,t) obeys to: 

= -Dk^ P{k, t), P{k, < = 0) = 1, (4.18) 
with the obvious solution P{k, t) — e^^^^*-, or, finally. 

This shows a purely diffusive behaviour: the probability to find the particle 
within a fixed volume centered at some point x goes smoothly to zero as i — ^ oo 
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for any x (runaway solution). The correlations of x reflect this behaviour too; 
for instance: 



r^{t) = {{x-xaf){t) = / (fix{x~ xof P{x~ XQ,t) = 6Dt, (4.20) 



showing that, on the average, the particle gets further and further away from 
the original point xq, but along a non-differentiable trajectory: f{t) cx -y/t, so 
the average velocity v — f{At)/At has no well-defined limit when At 0. 
This situation may change, however, if the motion of the particle is biased by 
an external force. Assume this force to be derived from a potential: F' = 
—dV/dx^. Then one can check that the tim e-inde pendent distribution Po{x) ~ 



exp[—(3V{x)] is a stationary solution to eq. ( 4.15 ) provided (3D = 1. Of course, 
this solution is acceptable as a probability density only if it is normalizable, 
which puts some constraints on the form of the potential. But assuming this to 
be the case, then Po{x) ^ e~^'^ represents an equilibrium distribution which is 
(asymptotically) reached by the system at large times |5g|l . O nce this is done. 



all the correlations become independent of time (unlike ( [4.20 )). This solution 



is a "fixed point" in the functional space of all (acceptable) distributions. 



Returning to our RGE ( 3.49 ), it should be clear by now that this is a functional 
Fokker-Planck equation which describes a random walk in the functional space 
of the colour fields a'^{x±). In this equation, t] plays the role of the "diffusion 
coefficient" , while v is like a "force term" , although this identification is somehow 
ambiguous since rj is itself a functional of a, so its derivatives can generate other 



contributions to the force term, as shown in eq. (4.7). (In the analogous problem 
of the Brownian motion, this would correspond to a diffusion coefficient which 
depends on x and has a tensorial structure: D — > Dij{x). This situation occurs, 
e.g., in the description of a random walk on a curved manifold [p9[.) In fact, 
it is more correct to identify the combination ^{Srj/Sar) ~ v a.s the effective 



"force term", since the remaining second-order differential operator in eq. (4.7) 
— which describes diffusion — is then Hermitian and positive semi-definite. 



A fixed point of the quantum evolution would be a solution W[a\ to eq. (3.49) 
which is normalizable and independent of "time" r. If such a solution existed, 
then the high energy limit of QCD scattering would be trivial (at least, within 
the present approximations): At sufficiently high energies, all the cross sections 



would become independent of energy (recall that r ~ Ins). The relation (4.8) 



between the coefficients in the RGE guarantees, however, that such a "fixed 



point" does not exist: The effec tive f orce in eq. (4/7) vanishes, and the corre- 



sponding evolution Hamiltonian (4.13) is just a kinetic operator, which describes 
pure diffusion. Wc thus expect gluon correlations to keep growing with r ~ In s 
even at asympt otical ly large energies. In Sect. I, we shall find approximate 



solutions to eq. ( 4.1 2| ) which show indeed such a behaviour 
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4.3 The Balitsky-Kovchegov equation 

If {0[a\ )t is any observable which can be computed as an average over a : 



{0[a])r 



V[a]0[a]Wr[al 



(4.21) 



(cf. eq. (2.41)), then its evolution with r is governed by the following equation: 



d 



(OH 



VaO[a 



1 



dWr[a] 
dr 



2 <5a?(a;_L) '5a^(2/±) 



0[c 



(4.22) 



where, in writing the second line, we have used eq. (4.12) for dWr/dr and then 
integrated twice by parts within the functional integral over a. 
Let us apply this to the 2-point function (2.26) of the Wilson lines in the funda- 
mental representation. We recall that, physically, this is the S'-matrix element 
for dipole-hadron scattering (cf. Sect. 2.2.2| ). A straightforward calculation 
yields (see for details): 

^2 



d_ 



{iv{V^,Vy))r = 



d zj_- 



{xj_ - yj_y 



(4.23) 



27r2 7 (a;_L - z_L)2(yj_ - z_l)2 

X (NMv^Vy) - tTiv:^v.Mv,%))^ 

This is the equation originally obtained by Balitsky |44|, within a quite different 
formalism : by an analysis of the quantum evolution of the dipole itself. 
Not that the above equation is not closed: It relates the 2-point function to 
the 4-point function {ti{V^V^)ti{V^Vy)) . One can similarly derive an evolution 
equation for the latter |44| , but this will in turn couple the 4-point function to 
a 6-point function, and so on. That is, eq. ( 4.23| ) is just the first in an infinite 
hierarchy of coupled equations [^ . 

A closed equati on can still be obtained in the large limit, in which the 4-point 
function in eq. (1^ factor izes: 

(tTiV^V,) triV^Vy))^ ^ (tr(T4^V;))^ HV^Vy))^ for TV, ^ oo. 
Then eq. ( [4.23| ) reduces to a closed equation for St{x±, y±) ~ (tr(Vjl/^))^/iVc : 



5 / X OLsNc 

— 6^(a;_L,2/_L) = 



(a;_L - y\_f 



27r^ J (a;_L — z_L)^(y_L — z_l)2 

X {S'r(xj^, y_L) - S'r(x_L, z_L)S'r(z_L, y^)} • (4.24) 

The same equation has been independently obtained by Kovchegov p5| within 
Mueller's dipole model (See also Ref. |Q for another derivation.) 



An important observation refers to the transverse kernel in eqs. (4.23) or (4.24): 
This is not the same as the original kernel /C(a;_L, y_L, zj_), eq. (4.11), of the RGE. 
Rather, this has been generated as 



/C(a;j_, a:_L, z_l) + /C(y_L, y_L, 2:_l) - 2/C(a;_L, y_L, z_l) = 



(a;_L - y\)" 



(x_L - z_l)2(?/j_ - zaJ^ 
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and has the remarkable feature to show a better infrared behaviour than eq. (4.11): 
When zj, a;_L , , the kernel above decreases like {x± — y±y/z^, so its integral 
over z± is actually finite. 

There is currently a large interest in the solutions to eq. (4.24), and significant 
progress has been achieved by combining analytic and numerical methods p9l 
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55| , |58| . The conclusions reached in this way are equivalent to those 
obtained from direct investigations of the RGE (3.49) ||l^, |5^ that we shall 
review in what follows. 



Approximate solutions to 

the Renormalization Group Equation 



We shall now construct approximate solutions to the RGE ( [4.12 ) and study 
their physical implications p7| . 

5.1 The mean field approximation 

As compared to the standard diffusion equation ( 4.15| ), the main complication 



with the RGE ( 4.12 ) com es fro m the fact that its kernel rj is itself dependent 
on a. In this respect, eq. ( [1.12| ) is similar to the following diffusion equation: 

dP(x,t) d „ , . d , . , , 

in which the diffusivities Dij (x) are allowed to depend upon the position x of 



the particle. This dependence makes eq. (5J) difficult to solve in general (i.e., 
for some arbitrary tensor field Dij{x)). But since a; is a random variable, with 
probability density P{x, t), a reasonable approximation is obtained by replacing 
Dij{x) in eq. (|5.lD by its expectation value: 



D^J{x) —> {D.,,{x))(t) = / £xP{x,t)D,j{x) = S,,D{t), (5.2) 



*j \ / \ u V //V / — J \ ' / \ ) — u 

which is independent of a:, but a function of time. We denote with a bar quanti- 
ties evaluated in this "mean field approximation" (MFA). In particular, P{x,t) 
is itself related to D{t), as the solution to the following approximate equation: 

dP(x,t) _ ^ 92 _ 

Thus, eq. ( ^.2[ ) is actually a self-consistent equation for D{t). Being homoge- 
neous in X, eq. ( ^.3[ ) is easily solved by Fourier transform, as in eqs. ( 1.18 )-( pl.l9| ). 
For the initial condition P{x,t = 0) = S'^^\x), one thus obtains: 
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By inserting this solution in eq. (5.2), one can compute the average there (as 
a functional of D{t)), and then solve the self-consistent equation for D{t), thus 
completely specifying the approximate solution (5.4). 

This is the strategy that we shall use to obtain approximate solutions to the 
functional diffusion equation ( 4.12| ). The corresponding MFA reads: 



dWr[a] 
dr 



with 



(5.5) 



(5.6) 



where the trivial colour structure in the l.h.s. follows from gauge symmetry. By 
the same argument, {v°'{xa_))t = 0, which is indeed consistent with the MFA 
[b.t) for 77 and the condition (^^). 

Eq. (5.5) is homogeneous in the functional variable ay{x±) (since its kernel 7r is 
independent of a), so it can be solved by functional Fourier analysis. This is the 
straightforward extension of the corresponding analysis for ordinary functions, 
and can be more rigourously introduced by using a discretized version of the 
3-dimensional configuration space (y,2:j^), as in eqs. (2.5) or (4.4). We write: 

5[a] = J D[n]e-'-f'^y-f'^^''^''y^'-'^'X^''^\ 

Wr[a] = /D[7r] e-^'''^y^'*'^^'^?(^^)""(^^)iy^[7r]. (5.7) 



By inserting this representation for Wr[a\ in eq. (5.5), and using 
6 



5a'^{x±) 



dy d zj_Tr'y{zj_)a'y{zj_) = 7rJ^(a;j_), 



one obtains the following equation for WViTi"] (compare to eq. ( 4.18 )): 
dWrM 



dr 



(5.8) 



(5.9) 



with the immediate solution (transverse coordinates are omitted, for simplicity): 

Wr [tt] = e" 5 lo dyiy-;^ Wo [k] . (5.10) 

The argument TTy of the initial weight function Wo[7r] has support only at y < 0. 
After insertion in eq. (5.7), this yields: 



Wr[oi] - / DI-k] e-'S-.^'^y<'^y e-h!odyiyKK = 5[a>]Wr[a^ 



(5.11) 
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with 6[a-^] defined in eq. (4.4) (this has been generated by the functional integral 
over TTy with y > t) and 



W,[a<] = AA.exp -- / dy 



Wo [a] 



(5.12) 



In this equation, yVo[Q^] is the original weight function at t = 0, and is a func- 
tional of the fi eld c ty with y < 0. (A/V is an irrelevant normalization facto r.) 
The solution (pll)-(p^) has the general structure anticipated in eqs. ( [4.2[) - 
(^!^). If the initial conditions are described by t he M V model, or an y other 
MFA, then Wo [a] is a Gaussian too (see, e.g., eq. ( 2.52| )), and eq. ( 5.12 ) can be 
rewritten as: 



Wr[a<] = TV^exp 



dy 



(5.13) 



For y < 0, the width 7y is specified by the initial conditions, while at positive 
rapidities < y < r, it is determi ned b y the quantum evolution, as we shall see. 
The fact that the weight function ( 5.13 ) is a Gaussian does not necessarily mean 
that the present approximation describes a system of independent colour sources 
(like the MV model). It just means that, in the MFA, all the correlations are 
encoded in the width of the Gaussian, or, equivalently, in the 2-point function 



(a°(a;_L)a^,(y_L))7 



(5.14) 



But this 2-point function contains also information on the higher-point cor- 
relations, although just in an averaged way, because it is determined by the 
following, non-linear, self-consistency equation: 



7r(a;±,y_L) = 



(2^ 



IC{x±,y±,zj_) 



(5.15) 



X (l + Sr{x±,y±) - Srix±,zj_) - Sriz±,y±)^, 



which follows from eqs. ( p.q ) and (3.51) together with the fact that, for a Gaus- 
sian weight function^^, 



sab 



m - 



-iTY{v\xi_)V{yA_)))^ = 6-'Sr{x^,y^), (5.16) 



with St a (non- linear) functional of 7y, to be constructed shortly. 
The correlation function (5.14) is local in y : colour sources located at dif- 
ferent space-time rapidities appear to be statistically independent. This is, of 
course, just an artifact of the MFA. The complete RGB generates correlations 
in rapidity, via the Wilson lines in its coefficients. But the only trace of th ese 
correlations in the MFA is the fact that the self-consistency equation ( 5.15| ) is 
non-local in y. 



Note that, as compared to cq. (2.26), St is now written in the adjoint representation. 
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To perform the average in eq. ( ^.16 ), we first derive an evolution equation for 
Sr, by using the corresponding equation (bTs) for Wr : 



d f 

-^Sr{x±,y±) = / D[a] 



D[a]Wr[a] 



(5.17) 



\lrix±,xj_) +7r(y±,y±) - 27r(a;_L,2;±)JS'r(a;_L, j/_l). 
(The functional derivatives of the Wilson lines have been evaluated as 

5a-{u^)Sa-{v^) ^'''^^^^^ = -9'NcTr{V^Vy)iS,, ~ (5,„)((5,. - Sy^), (5.18) 



where we have used eq. ( 3.53 ) and T'^T'^ — Nc-) Eq. ( 5.17 ) can be trivially inte- 
grated. To simplify the calculations, we assume homogeneity in the transverse 
plane within the hadron disk of radius R; then ^y{xj_,yj_) ~ "/y{xj_ — yj_) and 

S,{r^) = e-f'^-''o^^y[^>-("-)-^-('^-)l5o(r^) = e-f'^=K-(°-)-«-('--)l , (5.19) 
where r± = x± — y±, 

Ur±) = Ur±) + r dyjyir^), (5.20) 



- p-9'JVcKo(0_L)-«o(ri)] 



and in writing the second equality in (5.1£) we have assumed that the initial 
condition So{rj_) can be written in the form So{rj_] 



This is indeed the case for the weight function in eq. (5.13) — in particular, for 
the MV model, cf. eq. ( 2.6C| ) — , which yields : 



(,(){r±) = / dy7y(rj.). 



(5.21) 



By combining eqs. ( 5.15| ), (5.19) and (5.2C), one can finally rewrite the self- 
consistency equation as an evolution equation for ^t('"j_) : 



d^r{xj_ - yj_) 
dr 



1 f(Pz_i_ 



IC{xj_,y±,z±) 



(5.22) 



ttJ (27r)2 

X (l + Srix_i_ - y_i_) - Sr{x± - z±_) ~ Sr{zi_ ~ y±)j , 

with Sr{r±) = c^9'i^a[iT{o±)-irir±)] ^ anticipated, this equation is highly non 
linear in ^r- It is furthermore non-local in the transverse coordinates, but local 
in the "evolution time" r. (The original non- locality of eq. ( 5.15| ) in y has been 
now absorbed in the relation ( 5.20| ) between ^t- and 7y.) 

In the next sections, we shall develop further approximations, which rely on the 
kinematics and allow us to make progress with eq. ( ^.22 ). 
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5.2 Saturation scale and Kinematical Approximations 

Both the non-local and the non-linear structure of the evolution equation ( [5.22 ) 
depend crucially upon the behaviour of Srir^) with the transverse separation 
r±. From its definition (5.16), it is clear that >S'.r(^-L) ^ 1 as r_L — > for any 
T. Moreover, since a large dipole is strongly absorbed by a hadronic target, 
we expect that Sr{r±) 1 for sufficiently large r±, where what we mean by 
"sufficiently large" will generally depend on r. For instance, we have seen in 
Sect. |2j, within the MV model, that S'^(rj_) < 1 for rj_ > l/ Qs, w ith the 
saturation scale for gluons in the hadron wavefunction (cf. eq. ( ^.76| )). In that 
classical model, Qs was independent of energy, but in general we expect it to 
increase with r, because of the quantum evolution (cf. the discussion in Sect. 
1.4 and Sect. 5.3 below). At a formal level, this intimate connection between 
the strong absorbtion limit for a colour dipole and gluon saturation is based on 
the fact that, in both problems, the non-linear effects are encoded in Wilson 
lines. So, let us introduce the correlation length 1/Qs{t) of Sr{r±) : 



1, for rj_ < 1/Qs{t) 
0, for ri»l/g.(T) 



(5.23) 



which, as its notation suggests, will play also the role of the saturation scale. 
This behaviour of S't(?'_l), with an unique separation scale between a short-range 
regime and a long-range one, is confirmed by numerical studies of the Kovchegov 
equation, which also show a rapid increase of Qs with r ^ ^ p8[ . 
Eq. (5.23), together with the expression (5.1£) for 5r(r_L) in the MFA, imply 
the following condition: 



'NcMO±)-^Ar±)] ^ I for r^^l/Q.ir), 



(5.24) 



that we shall use later to obtain an estimate for Qs{t). 

An external probe with transverse momentum k± will measure correlations in 
the hadron over a typical transverse size rj_ ~ l/kj_. Thus, short distances r± <gC 
1/Qs{t) correspond to high momenta, k± 3> Qs(r), while large separations 
rj_ ^ 1/Qs{t) correspond to low momenta k± <C Qs(r). In what follows, we 
shall not aim at a precise description of the physics around the saturation scale, 
but rather focus on the two limiting regimes — high-Zcj^ and low-fc_L — and 
perform appropriate simplifications on the evolution equation ( 5.22| ). 



a) High-k± . It is convenient to rewrite eq. ( 5.19|) in momentum space as: 



(27r)2 



(5.25) 



For rj_ <C 1/Qs{t), the integral over p± is dominated by momenta within the 
range Qs (r) <Cp_L <Cl/rj^ . This holds to leading transvers^log accuracy: In this 
range, ^r{p±) ^ ^Ip\ (up to logs), so the integral over p_L produces the large 
logarithm \T^(\ j r\Q\(Ty) . To the same logarithmic accuracy, one can expand 



61 



the exponential in ( 5.25| ) in powers of p± ■ rj_, like in eq. (2.64), and thus obtain: 

l/r5 



exp<! 



(27r)2 



(5.26) 



When extrapolated to r± ^ 1/Qs{t), this expression gives us the following 
estimate for the correlation length 1/Qs{t) (cf. eq. (5.24|)) : 



a,N, 



j dp\p\^r{Pl^)- 



(5.27) 



For ri_ <C 1/Qs{t), where eq. ( 5.26 ) is strictly valid, the exponential there can 
be expanded to lowest order: 



Sr{ri_) 



1 



c 2 

r 



(5.28) 



(The ultraviolet cutoff l/r_L is implicit in the momentum representation of 
^r(O).) By inserting this into ( 5.22 ), we obtain a linear evolution equation 
for ^r(r_L) : 



dS,T{xi_ - yi_) 
dr 



tPz_i (x* — z*)(y' — z') 



(5.29) 



(2^)2 (a;i-zi)2(y^-z^)2 
X [{xi_ - vi_f - ixi_ - z^f - [z^ - 2/±)') Vier(O) 



Thus, the short-distance approximation is automatically a linear, or weak-field, 
approximation. This is to be expected since, at high fcj^, the gluon density is 
low. 

To perform the integral over z± in eq. (5.2S), it is useful to recall eq. (3.48) and 
then notice that, within the integrand of (|5.29 ), one can effectively replace: 



1 (x* - z*)(y* - 
(2^ {x^-z^Y{y^-z^Y 



(The additional terms in the r.h.s. are (5- functions at zj^ — x±_ or zj^ — 
which vanish when multiplied by the remaining factor in (5.29).) By using this, 
together with a couple of integrations by parts w.r.t. z±, and a Fourier transform 
to momentum space, we finally obtain the following evolution equation: 



dr 



a,N, 



dpi , . 



(5.30) 



for the quantity: 



(5.31) 
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which, physically, is the 2-point function of the colour charge density in the 
transverse plane p°'[x±_) : 

{p''{x^)p\y^))r = S'^'fir{xi_ - y^), p'^ix^) = [ dyp^(.T^). (5.32) 



The initial condition for eq. (5.3C) can be taken from the MV model: pT^{kj_) = 
jjLA for T = 0, cf. eq. (2.51). This initial condition is independent of k±^ and, 
together with eq. ( ^.30 ), it implies that iiT{kj_) remains a rather slowly varying 
function of k± in this high momentum regime. This will be manifest on the 
solutions to eq. ( |5.30| ) that we shall write in the next subsection. 

b) Low-k± . For large distances r± 1/Qs{t), Sr{r±) ^ 1, and the 2-point 
functions of the Wilson lines can be simply neglected in the self-consistency 
equations ( |5T5| ) or ( |5^ ) 0, H]. Eq. ( |5^ ) then simphfies to (see also 
eq. dH)) 



d^T{x± - yj_) 
dr 



^ Id'z^diix^l 



\z±)dl{yi 



-Vi 



\y±) 



or in momentum space (cf. eq. (5.20)) 



lr{kl_) = 



1 1 



dr 



(5.33) 



(5.34) 



This is not an equation anylonger, but rather an explicit, and rather simple, ex- 
pression for the propagator 7r (A:j^) of the fields a : this is just the 2-dimensional 
Coulomb propagator. 

Rem arkab ly, the QCD coupling constant g has dropped out from eqs 



(5.33) 



and ( ^.34). ( This should be contrasted with the corresponding equation at high 
fc_L, eq. ( |S.29 ), whose r.h.s. is exphcitly proportional to as = g'^ /^-k.) The same 
property holds then for the corresponding mean-field Hamiltonian (cf. eq. (5.5)): 



low — /C_L 



1 

2^ 



(Pki 
(M 
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kl 



5a%{ki_)5a%{-kx_y 



(5.35) 



which is quite remarkable since at low k± we are effectively in a strong coupling 
reg ime (in the sense that the COV-gauge fields are strong: a" ^ 1/g ; see Sect. 
5.4). If g nevertheless drops out in this limit, it is because of the special way it 
enters the evolution Hamiltonian: via the exponent of the Wilson lines. That 
is, the relevant degrees of freedom in the non-linear regime are not the (strong) 
colour fields by themselves, but rather the Wilson lines built with these fields. 
The Wilson lines are rapidly oscillating over distances r±_ ^ 1/Qs{t) (since 
their exponent is of order one, and the typical scale for variations is 1/Qs(r)), 
and thus average to zero ("random phase approximation"). 
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For what follows, it is useful to summarize the previous kinematical approxima- 



tions into the following, factorized, form for the weight function (5.13), which 
is most conveniently written as a weight function foij^ Py{kj_) = kj^a^{k±) : 

. r Qs{y) . 

,.;lowr 1 J f ^ f rf"fc± pg(fc^)pg(-fc±) I 

[p] ^ AA.expj-- J dy J ^ 1, (5.37) 



^ -°° Q4y) ^ 

In writing this equation, we have separated, for each rapidity y, the low-momentum 
(fcj^ < Qs{y)) modes of p fro m th e high-momentum (k± > Qs{y)) ones, we 



have used the approximation ( 5.34 ) for the width of the Gaussian at low mo- 



men ta, an d we have written Ay(fcj^) = dpy{kj^)/dy, with prikj_) determined by 



eq. ( 5.30 ), at high momenta. Note that the modes with k± ~ Qs(y) are not 
correctly described by the present approximations, but we shall assume that 
they give only small contributions to the quantities to be computed below. 

5.3 High— A;^ : Recovering the perturbative evolution 

We now consider the implications of eqs. ( |5.30| ) and ( |5.38|) for the physics at high 
t ransv erse momenta k± ^ Qs{t). To this aim, we compute the gluon density 
( ^.23 ) in this low density regime, where one can use the linear approximation 



jF+^(fc) ~ {ikP/Ji\_) p{k) . The calculation is similar to that already performed in 
eqs. (pt5^)-(^^ Specifically, by using (cf. eq. (|3|)) : 



{p^^{x^)p\,{y^))r^5'^'5{y-Y')\y{x^-V^). Ay(r^) = (5.39) 



dy 



one eventually obtains: 



xG(x,g2) ^ ^^^R^ I ^Prik^)- (5.41) 

Qi(r} ^ 

Note the lower limit Qs{t) in the integral giving xG(x, Q^) : for ^ Q1{t), 
and to leading transverse-log accuracy, it is sufficient to consider the contribu- 
tion of the high-fcj^ modes of p to the gluon distribution. We shall check later 
that the corresponding contribution of the modes with fcj^ <C Qs{t) is infrared 

^''This is the colour charge density in the COV-gauge, but we omit the tilde symbol on p, 
to simplify writing. 
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finite, although subleading as compared to eq. (5.41) IT^, |5^. This cures the 



infrared problem that we have faced in the classical calculation of Sects. 2.4-2.5 



Physically, HT{k±_) plays the same role as ha in the MV model: It measures 
the density of the colour sources in the transverse plane, and, in the linear 
regime at high-fc_L, it is also proportional to the unintegrated gluon distribu- 
tion: fi-riQ^) oc dxG{x,Q'^)/d\nQ^. But unlike ha, which is constant for a 
given atomic number A, iJ,Tik±) has non-trivial dependences u pon b oth r and 
k±, as determined by its quantum evolution according to eq. ( 5.30| ). The de- 
pendence on T describes the increase in the density of the colour sources via 
soft gluon radiation. The dependence on k± corresponds in coordinate space to 
correlations in the transverse plane, which occur via the exchange of quantum 
gluo ns (se e Fig. [T7|.a). 

Eq. (Is^Ol ) can be recognized as the standard, linear evolution equation in the 
double- logarithmic approximation (DLA) ||], i.e., in the limit in which BFKL 
and DGLAP coincide with each other. (In this limit, o nly th e first, "real", term 
must be retained in the r.h.s. of the BFKL equation (3.59); for k± 3> p±, this 
term reduces indeed to that in eq. (pel).) The emergence of DLA is natural, 
given the approximations performed in deriving eq. (5.30): we have kept only 
terms of leading-log accuracy in both t = ln(l/x) (in the construction of the 
effective theory), and ln(/c^/Q^(r)) (in the short-range expansion at high k±). 
Eqs. (|3^) and (|1^) imply the more standard form of the DLA equation B : 



92 



drdlnQ^ 



xG(x, g2 



a,N, 



xG(x, Q2). 



(5.42) 



At large r and/or Q^, the solution to this equation increases like (with ag 
UsNc/tt and Qq some scale of reference) [|j 



xG(x,Q^) oc exp<^ 2Ja,T ln(Q2/(3o) 



(5.43) 



where we have assumed as to be independent of 

o account, by writing a. 
then the dependence of the solution upon gets softer Q 



running of the coupling into account, by writing as{Q^) — bo/ \n{Q'^ /Aq^jj), 



If instead one takes the 

QCDJ 



xG(x, Q^) oc exp 




r In 



ln(Q2/A2 



QCDJ 



(5.44) 



In any case, eqs. ( ^.43 ) and ( 5.44 ) show that, at high transverse momenta 



Q ^ Qsi''')i the gluon distribution xG(x, Q^) grows rapidly with r. This 
is the standard picture of parton evolution, which, if extrapolated to arbitrar- 
ily high energies, would predict violations of the unitarity boundf^ But from 
the previous analysis, we know that the approximations leading to eq. (5.42) 
will break down at sufficiently large energies, where the non-linear effects in 



^*Note that altli niigli slnw pr than for the BFKL solution (3.6), the growth with t of the 
DLA solution (5.43) or (p.44[) is still faster than that of any power of r ~ Ins. 
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the quantum evolution cannot be neglected anylonger. Alternatively, for fixed 
rapidity r, the linear approximation breaks down at low transverse momenta 
k± <C Qs{t), with Qs{t) the saturation scale. An estimate for this scale has 
been given in eq. ( [S.27 ), which, together with eqs. (5.31) and (5.41), implies: 



^-Mr(p_L) = 
PI 



na.Nr, 1 



-jxG(x,Qf(r)). (5.45) 



1 i?2 



By further combining this result with eq. (5.43) or (5.44), one can deduce the 
T-dependence of the saturation scale in the DLA. One thus obtains: 



(fixed coupling), 



(5.46) 



_ a2 ^V2boTlnT 



(running coupling) . 



(5.47) 



and, respectively, 

Eq. dsl^ ) (or ( p7| )) defines a curve in the t — k± plane, which divides this 
plane in two (see Fig. |l8|) : Points on its right are effectively in the high 
momentum regime; they correspond to a dilute gas of weakly correlated colour 
sources whose density is rapidly increasing with r. Points on the left of the 
saturation curve correspond to the low momentum regime, to be discussed in 
the next subsection. 



5.4 Low— A;^ : Coulomb gas and gluon saturation 

We finally turn to the most interesting physical regime, that of the non-linear 
physics at small transverse momenta k± <^ Qs{t) (with k± 3> Aqcd, though), 
whose understanding was a main motivation for all the previous developments. 
Within the effective theory, the low-momentum modes of the colour source 



are described by the weight function W^™ , eq. ( 5.37 ), which is equivalently 
rewritten as (cf. Fig. |8|): 

with f{k±) — the rapidity at which the saturation momentum is equal to k± : 

Qlifik^)) = kl. (5.49) 



There are several noteworthy features about eq. (^ 
i) This describes a Coulomb gas, i.e., a system of colour charges interacting via 
long-range Coulomb forces. The colour source Py{x±) at x± feels the Coulomb 
field ay{x±) created at x± by all the other sources: 



p^{k^)p'^{-k^ 



,2 I ryV--L;\--L| 



p''{xi_){x^\ —2- |y_L)Py(yj.) = / p'^{x±)a^ixj^ 
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T(kJ 



Linear Evolution 
(Dilute Gas) 




DGLAP 



Figure 18: A "phase-diagram" of the various regions for evolution in the t — fcj_ 
plane. 



The fact that the charge-charge correlator appears to vanish when fc_L ^ is 
in agreement with gauge symmetry: The colour source Py{k±) at low k± is an 
induced source, whose global strength must vanish: 



x±,y±'Jy,y' 



= 0. (5.50) 



ii) The colour charge correlations are local in rapidity : the Coulomb forces 
couple only sources located in the same layer of y (or x^). At low-Zcj^, this 
property is not just an artifact of the MFA, but rather has a deep physical 
meaning: In the quantum evolution, the colour sources at different rapidities 
get correlated with each other because of the presence of Wilson lines in the 
evolution Hamiltonian (4.13). But these correlations are washed out on a large 
scale r± ^ 1/Qs{t), on which the Wilson lines average to zero. In particular, 
this explains why the width cx A:^ of the Gaussian (5.48) is independent of 
the initial conditions at r ~ 0. (By contrast, at high momenta, the width 

) is sensitive to the initial conditions, since 



dpy{k±)/d'y in eq. (5.3 



determined by solving eq. ( 5.30[ ).) 
iii) According to eq. ( 5.48| ), the low-momentum modes of p are uniformly dis- 
tributed in rapidity, within the interval f{k±) < y < r. It follows that the 
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integrated quantity: 

/^r(fc_L) = 



/r2 



' TT 



(5.51) 



which measures the density of sources (with given k±) in the transverse plane, 
grows only linearly with r, that is, logarithmically with the energy. This is 
to be contrasted with the strong, quas i-expo nenti al, increase of /iT-(fcj^) in the 
high-momentum regime (cf. eqs. ( 5.45 ) and ( 5.44|) ). We conclude that, at low 
momenta k± <^ Qs{t), the colour sources saturate, because of the strong non- 
linear effects in the quantum evolution. 

iv) The saturated sources form the outermost layers of the hadron in the longitu- 
dinal direction: for given kj_, they are located at > XqC'^^'^^K In particular, 



(5.52) 



is the longitudinal extent of the saturated part of the hadron, in units of rapidity 
(for modes with tranverse momentum k±). In writing ( 5.52| ), we have used the 
DLA estimate ( 5.46 ) for the r-dependence of the saturation scale, 
v) Note the factor l/ag in the r.h.s. of ( 5.52 ); this implies that, at saturation, the 
integrated charge density p°'{x±) has typically large amplitudes: p ^ \/{pp} ~ 
1/g. The same is therefore true for the COV-gauge field a"'{x±) : a ^ 1/g. 

Since the colour sources at low-fc_L are saturated, there should be no surprise 
that the gluons emitted by these sources are saturated as well, and this indepen- 
dently of their mutual interactions (i.e., of the non-linear effects in the classical 
Yang-Mills equations). Indeed, a quasi- Abelian calculation of the gluon distribu- 
tion, based on the linearized solution JF+^(fc) « {ik^ /k']_)p, yields the following 
gluon density (cf. eqs. ( |5.40| ) and eq. ( 5.51 )) : 



1 1 



167r*c a< 



In- 



kl 



(5.53) 



which already exhibits saturation ! In fact, as argued in Refs. |1^, the only 
effect of the non-linearities in the classical Yang-Mills equations in this low- 
fc_L regime is to modify the overall normalization of the linear-order result. In 
anticipation of this, we have inserted in eq. ( ^.55| ) a corrective factor c, which 
cannot be accurately determined in the present approximations (since sensitive 
to the physics around Qs), but should be smaller than one (although not much 
smaller) . 

Note the striking similarity between eq. ( ^.53 ) and the corresponding predic- 
tion (2.71) of the classical MV model. Despite of the differences in the physical 
mechanism leading to saturation — non-linear quantum evolution for eq. ( 5. 53] ), 
as opposed to non-linear classical dynamics for eq. (2.71) — , the final results 



look very much the same. So, the earlier discussion of cq. (2.71) can be imme- 
diately adapted to eq. (5.53), after replacing A —i- s : Eq. ( ^.53 ) shows marginal 
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A 



QCD 



Figure 19: The gluon phase-density Mrik^) in the effective theory plotted as a 
function of kj_ for two values of r = ln(l/x). 



saturation (in the sense of a logarithmic increase only) with both s and l/k'j_, 
with a typical amplitude of order l/ag. This is illustrated in Fig. [l9| , which 
should be compared to Fig. ^ (The high-fcj^ behaviour in Fig. |19| is taken 
from eq. ( ^.40D .) 



Aside from saturation, eq. (5.53) has also other important consequences, which 
all reflect the proportionality to the rapidity window[^ t — f(fcj_), eq. (5.52): 
a) Scaling. The gluon density at saturation depends upon the energy s and the 
transverse momentum k± only via the scaling variable 



T: 



\{r)/kl 



(5.54) 



A similar scaling is observed in the solutions to the Kovchegov equation 
p6| , |58{ . As mentioned in the Introduction, such a scaling has been actually 
observed in DIS at HERA Q. 

h) Universality. Eq. ( 5.53| ) is only weakly sensitive — via its logarithmic de- 
pendence upon the saturation scale — to the initial conditions for quantum 
evolution, and therefore to the specific properties of t he ha dron under consid- 
eration (e.g., its size and atomic number). Thus, eq. ( 5.53| ) not only provides 
arguments in the favour of hadron universality at high energy, but also predicts 
what should be the pattern of its violation. 

^^These properties are therefore generic: They hold for any quantity which receives his 
dominant contributions from the saturated gluons. 
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The gluon distribution xG(x, Q^) at ^ Qs(''') is immediately obtained by 
integration in eq. (5.53): 



xG(x, Q' 



N^-l 1 



kl 



■^R'Q^ [HQ^,{t)/Q^) + 1] 



(5.55) 



Note that, since ijLr{k±) ^ k\ in the saturation regime (cf. eq. ( 5.51 )), the 
above integral is almost insensitive to the soft modes kj_ ^ ^qcd- This has 
allowed us to extend the integration down to fcj^ = without loss of accuracy. 
As anticipated, the phenomenon of saturation reduces the sensitivity of physical 
quantities to the infrared gauge fields, thus making the weak coupling expansion 
reliable. (In Ref. a similar conclusion is drawn on the basis of Kovchegov 
equation.) If extrapolated up to Q Q^, eq. (5.55) yields 



xG(x,Q2(^)) 



1 1 



167r^ 



(5.56) 



in rough agreeme nt wi th the corres ponding extrapolation from the high momen- 
tum regime, eq. ( 5.45 ). Eq. ( 5.56 ) gives also the contribution of the saturated 
modes to the gluon distribution at momenta Q > Qs{t). But for very high 
momenta, Q 3> Qs{t), the dominant contribution comes form the hard modes 
{Qs <^ k± ^ Q), and is given by eq. ( 5.41 ). 



As a final application, let us compute the 2-point ftmction S'r(?'_L) of the Wilson 
lines for large distances r± ^ 1/Qs{t). This is interesting for at least two 
reasons: It shows how the unitarity limit is reached for the scattering of a large 
colour dipole off the hadron, and it allows us to check a posteriori the consistency 
of the "random phase approximation" that we have used at low kj_ . 
To this aim, we rewrite eq. (5.25) as 



Srir^) ~ exp<^ - 



s(y) 



dy 



(5.57) 



where we have anticipated that the main contribution comes from the saturated 
modes, for which = ^/i'^P±)j cf. eq. ( |5.34 ). The integral over p± is now 
infrared finite (as opposed to the MV model: compare to eq. ( 2.64 )), and to 
leading log accuracy can be evaluated as: 



s{y} 



[ 


1 




J (2^)2 







(5.58) 



The result can be understood as follows: as long as l/r± » Qsiy), or y < T{r±), 
Qip±-r± ^ 2^ and the integral vanishes. Bur for y > f(rj_), or l/r± <C Qs{y), the 
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integrals corresponding to the two terms in the brackets are cut off at different 
ultraviolet scales: Qs{y) for the first term, and l/r± for the second one. Their 
difference gives the log in the r.h.s. By also using ln(Q^(y)r^) = 4:as{y — f{r±)), 
cf. eq. (5.52), and performing the integral over y, we finally deduce: 



Srir±) ~ exp{-2a2(r-f(ri))2} 



exp 



1 r 



HQl{r)rj 



(5.59) 



whic h coi ncides with the result obtained from the Kovchegov equation Q . 
Eq. (5.5!;) shows that the correlator of the Wilson lines is rapidly decreasing 



when Qj.{T)r^ S> 1, so that the RPA is indeed justified, at least as a mean field 
approximation. 

More details and further applications of the mean field approximation will be 
presented in Ref. where the results obtained in this way will be also 

compared to the corresponding predictions of the Kovchegov equation. It would 
be also interesting (especially in view of applications to phenomenology) to take 
into account the transverse inhomogeneity of the hadron (i.e., the dependence 
upon the impact parameter in the transverse plane). This can be done already 
in the framework of the MFA, but, more generally, it would be important to 
understand the limitations of the latter, and to be able to solve the complete 
RGE. This might be done, for instance, via numerical simulations on a lattice. 
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